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SECTION  I 


FAILURE  OF  SIMPLE  STRUCTURAL  ELEMENTS  SUBJECTED 
TO  IMPULSIVE  LOADS 


l. 1  Introduction 

Failure  of  a structural  element  may  be  defined  in  manv 
giving  a very  wide  range  of  choices  for  definition 
ffifpi  element  may  be  considered  as  hc.v5ng 

exceeded  some  load  or  stress  level,  some 
local  distortion  or  some  large  deflection  due  to  a buckline 
phenomena.  To  further  complicate  a failure  definition  various 
kinds  of  loads  may  produce  the  stress  level  or  distortion 
necessary  for  failure.  Failure  of  a structural  element  is 
Jui^P^Sn®  f term  and  a given  definition  must  be 

fSiSre  failure.  For  this  discussion 

^ ture  shall  be  defined  as  a condition  when  actual  rupture. 

racture,  or  breakage  occurs  at  some  point  in  the  material 
and  renders  the  structural  element  incapable  of  resisting  anv 
further  loading.  This  discussion  shall  be  x^estricted  to^fail- 

such  as  beams,  pla?esr^ods!  ete 

and  will  not  be  concerned  with  failure  of  a large  or  comnlicated 
structure  comprising  many  of  the  basic  elements?  Plicated 

ac,  4-vo  response  mechanism  of  a structural  element  may  be  defined 

deformation  leading  to  the 

associated  with  blast  loading  are  usually  one  of  two  types,  (1) 
a static  mechanism  where  the  initial  deflection  shape  is  con- 
tinued on  through  failure,  or  (2)  a traveling  plastL  hinge 

some^stIL”Sf°iha'Tf  ® deflected  shape  Continuously  for 

^ deformation  process.  The  response  mechanism 
^ structural  element  is  not  only  dependent  on  the 
general  properties  of  the  structure  but  is  also^s?rongirdepend- 

thp  h?pc=ih^  magnitude  of  the  peak  overpressure  and  duration  of 
the  blast  or  impulsive  loading.  “t-j-on  oi 

with  defined  in  a manner  associated 

anSm^^  Tf  Placed  on  the  structure  by  the  response  mech- 

in.teAsiL^  St\rL^SCrarSCniSiC\CuS  bC  CalL^rj  SLio^^"^ 
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In  order  that  a failure  load  may  be  predicted,  it  is 
necessary  to  know  both  the  mechanism  and  mode  of  failure. 

The  response  mechanism  determines  the  irreversible  path  taken 
by  the  structural  element  and  the  energy  required  or  work 
done  in  moving  the  element  through  that  path  may  be  predicted 
quite  readily.  Us,ing  the  failure  mode  based  on  knowledge  gained 
from  the  response  mechanism  a failure  criteria  based  on  a maxi- 
mum. stress  or  deformation  may  be  applied  for  that  particular 
failure  mode.  Since  complete  failure  of  structur’al  elements 
is  rather  complicated  it  is  best  that  a discussion  of  the  entire 
assumed  failure  process  and  resulting  analvsis  be  given  under 
one  general  heading  entitled  plates,  beams,  etc.  In  this  light 
the  remainder  of  this  section  will  be  devoted  to  a discussion 
of  plates  and  beams  subjected  to  impulsive  loads. 


1.2  Flat  Plates 

The  response  of  flat  plates  when  subjected  to  impulsive 
loads  resulting  from  blast  loadings,  as  mentioned  previously, 
is  usually  divided  into  two  types  of  mechanisms:  (1)  the  static 

case  where  the  plate  deforms  in  some  continuous  uniform  shape 
as  shown  in  Figure  la,  and  (2)  a plastic  hinge  case  where  the 
initial  deflection  is  by  propagation  of  stress  wave  or  plastic 
hinges  as  shown  in  Figure  lb.  For  plates  with  high  short  side 

Traveling  Plastic 
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dimension  to  thickness  ratios,  where  bending  is  negligible 
compared  to  the  membrane  forces,  the  resultant  failure  mode 
for^the  static  mechanism  is  tension.  The  location  of  the 
failure  is  dependent  on  the  plate  boundary  conditions  and 
distribution  of  load.  The  plastic  hinge  mechanism  operating 
in  the  same  plate  is _ characterized  by  a tensile  stress  wave 
(plastic  hinge)  originating  at  the  boundaries  and  moving  toward 
the  center  of  the  plate.  The  failure  mode  is  one  of  tension 
which  may  occur  near  the  edges  of  the  plate  in  the  early  s ^age 
of  plate  deformation  or  at  other  positions  in  the  plate  where 
high  strains  result  from  the  initial  transmission  and/or 
reflections  of  the  stress  wave.  Generally,  the  dividing  re- 
gion governing  which  mechanism  occurs  is  defined  in  terms  of 
the  characteristic  period  of  the  elastic  plate  frequency  and 
the  positive  phase  Of  the  pressure  pulse.  If  the  positive 
pressure  phase  of  the  pressure  pulse  is  greater  than  one- 
quarter  of  the  characteristic  period  of  the  plate,  then  the 
static  mechanism  is  assumed  to  operate  and,  if  the  positive 
pressure  phase  is  less  than  one— quarter  of  the  characteristic 
P^^iod , then  the  plastic  hinge  mechanism  is  the  governing 
mechanism.  However,  the  demarcation  is  not  distinct  and,  for 
the  region  where  a quarter  of  the  characterist5  c period  and 
the  positive  pressure  phase  are  approximately  equal,  either 
mechanism  may  occur  and  the  resulting  deformation  and  failure 
is  strongly  dependent  on  the  shape  of  the  pressure-time  curve. 
General^ discussion  of  the  response  of  circular  plates  subjected 
to  impulsive  loadings  is  given  by  Cristescu^l^  and  some  later 
work  IS  described  by  Johnson(2)  and  Abrahamson^ 3 ) . The  response 
of  square  or  rectangular  flat  plates  is  a much  more  complicated 
phenomena,  and  several  computer  studies  such  as  those  by  Kay(*+) 
and  Huffington  ) are  available.  Such  preliminary  results  on 
response  of  flat  plates  subject  to  fuel-air  explosive  blast 
waves  have  been  reported  by  the  authors^^^  and  a technical  re- 
port concerned  with  this  subject  is  in  progress;  therefore, 
further  aiscussi.on  on  response  of  flat  plates  will  be  omitted 
from  this  report. 


1.3  Previous  Analyses  of  Beams 

Seme  early  work  on  beam  response,  resulting  in  plastic  de- 
formation, IS  that  of  Duwez,  et  al^^^  and  is  based  on  trans- 
verse impact  of  a concentrated  load  on  an  infinitely  long  beam. 
An  elementary  analysis  of  a free  ended  beam  with  central  trans- 
verse impact  was  presented  by  Lee  and  Symonds(8)  at  about  the 
same  time  and  the  analysis  is  based  on  a rigid-perf ectly 
plastic  material.  This  assumption  makes  use  of  the  fact  that 
when  a large  amount  of  kinetic  energy  is  available  to  produce 
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plastic  deformation  the  plastic  work  in  the  beam  far  exceeds 
the  elastic  work.  This  assumption  used  in  References  7 and 
8 also  leads  to  the  use  of  a localized  plastic  deformation, 
commonly  called  a plastic  hinge,  as  a means  of  describing  the 
deflection  process.  The  plastic  hinge  mechanism  is  not  re- 
stricted to  dynamic  analyses  and  a general  discussion  of  the 
mechanism  and  several  applications  to  static  problems  is  given 
by  Timoshenko  and  Gere^^^.  Due  to  the  importance  of  the  assump- 
tions made  in  the  use  of  a plastic  hinge  mechanism  a general 
development  will  be  given. 

if  the  stress  in  a beam  remains  below  the  proportional 
limit,  the  beam  behavior  is  independent  of  the  shape  of  the 
plastic  portion  of  the  stress-strain  curve  of  the  material. 

When  a beam  is  loaded  in  such  a manner  that  the  stresses  ex- 
ceed the  proportional  limit  of  the  material  the  behavior  of 
the  beam  in  bending  is  then  termed  inelastic  or  plastic  and 
is  very  dependent  on  the  inelastic  portion  of  stress-strain 
curve  for  the  material. 

The  basic  assumption  used  in  this  discussion  is  that  once 
the  entire  portion  of  a bea.m  cross  section  exceeds  the  propor- 
tional limit  of  the  material  a plastic  hinge  forms  and  any 

increase  in  loading  will  cause  unrestricted  rotation 
of  the  beam  at  this  point.  The  plastic  moment  required  to  pro- 
duce^ the  plastic  hinge  is_dependent  on  the  shape  of  the  stress- 
strain  curve  of  the  material.  A discussion  of  the  various  stress- 
strain  curve  shapes  and  their  resultant  effect  on  the  plastic 
moment  and  hinge  will  not  be  given  here.  The  analysis  as  pre- 
sented in  this  report  will  incorporate  the  use  of  the  rigid- 
pi^stic  model  for  ductile  materials  as  shown  in  Figure  2.  Non- 
ductile  or  brittlematerials  tend  to  act  completely  elastic 
and  an  analysis  using  a plastic  hinge  is  not  applicable.  An 
elastic _ analysis  may  be  found  in  a discussion  of  classical  beam 
theory  in  texts  such  as  Reference  9. 

The  plastic  moment  for  a beam  is  also  dependent  on  the 
cross-sectional  shape  of  the  beam  and  is  readily  determined  in 
terms  of  the  yield  moment.  The  yield  moment  is  the  moment  re- 
quired to  produce  the  yield  stress  at  the  outermost  fiber  of 
the  beam  cross  section.  Using  the  stress  distribution  of 
Figure  3a  and  integrating  the  first  moment  of  stress  over  the 
beam  cross  section,  the  yield  moment  is  found  to  be 


(1) 


Strain 


Figure  2.  Idealized  Stress-Strain  Curve 

to  the  outermost  fiber.  uJing^F-gurraWhe'^SlaSr®'” 

may  then  be  determined  as  plastic  moment 


/ avydA 


Neutral 


Figure  3.  Bending  Stress  Distribut: 


plastic ^ moment  to  the  yield  moment  is  solely 
sh^p^fac^o?!  section  shape  and  is  defined  as  the  beam 

^ • (3) 

The  value  of  f for  various  shapes  is  given  below  in  Table  1. 

table  1.  beam  shape  factor 


Shape 

f 

2.00 

1.70 

1.50 

I 

1.15 

Using  Equations  (1)  and  (3)  the  plastic  moment  is  written 
Mp  = fay|  = fOyZ  , 

Where  Z is  the  section  modulus  defined  as 
Z = - . 

c (4) 

fo-PTTi  Stress  in  the  beam  is  such  that  plastic  hinees 

form,  th(=>  beam  is  assumed  to  act  as  a of  inges 

s2ppo?tJd°SLris^ioIded  stSticalif®;  fxamplet^/rsimple 

at  the  .idspan  and 
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the  midpoint  and  the  natural  hinges  at  the  ends  as  shown  in 

1*.  ^ beam  fixed  at  both  ends,  plastic  hinges  form 

at  the  fixed  ends  as  well  as  at  the  midpoint.  If  true  end 
fixity  occurs  such  that  the  end  points  do  not  move  axially, 
then  tension  is  induced  in  the  beam  due  to  the  deflection  and 
increased  length  of  the  beam. 


p(t) 


Hinges 


Plastic 

Hinge 


Figure  4.  Static  Beam  Failure  Due  to  Plastic  Hinge 

, The  dynamic  response  of  a beam  to  impulsive  loads  may  be 
described  in  terms  of  plastic  hinges  and  for  the  case  of  low 
overpressure  the  mechanism  may  be  the  same  as  the  static  analy- 
s.  However,  for  high  peak  overpressures  a traveling  plastic 
hinge , _ as  shown  in  Figure  5 , forms  and  moves  away  from  the  beam 
'4.  intersection  of  the  two  moving  plastic  hinges 
IS  dependent  on  the  beam  end  fixity.  Once  the  traveling  hinges 
intersect  then  the  response  is  the  same  as  that  for  the  static 
case.  The  sequence  for  this  mechanism  is  shown  in  Figure  5. 

condition  governing  which  mechanism  is  in 
IS  based  on  the  value  of  peak  overpressure  of  the  pressure 
time  function,  and  Abrahamson,  et  al^^^  show  that  if  the  peak 
overpressure  p < 3p^,  where  p^  is  the  static  colLpL  p^essSre, 
then  the  static  mechanism  is  valid.  Thus,  for  p > Sp^  a travel- 

^ development  and  discussion  of  equations 
motion  for  both  the  traveling  and  static  mechanism  in\  beam 


I ' ' 


neglecting  axial  tension  is  given  in  Reference  3.  A general 
discussion  of  the  results  of  this  reference  follows. 


Traveling  Plastic  Hinges 


I t=t. 


Fixed  Hinges 


Figure  5.  Moving  Plastic  Hinge  Mechanism 

The  analysis  of  beam  response  as  given  in  Reference  3 is 
divided  into  the  classical  static  and  plastic  hinge  mechanisms 
The  condition  governing  which  mechanism  holds  is  based  on  the 
peak  overpressure  of  the  loading  function  and  is  specified  in 
terms  of  the  static  collapse  pressure  pg.  The  static  collapse 

pressure  is  determined  from  a static  analysis  assuming  that  a 
uniform  pressure  exists  over  the  entire  span  2L  and  plastic 


hinges  are  operating  with  the  only  resisting  force  being  the 

plastic  moment  Mp.  Under  these  assumptions  the  static  collapse 
pressure 


Ps  = 4Mp/L2 


for  the  fixed  end  beam  and 


Ps  = 2Mp/L2  , 


for  the  simply  supported  beam.  Using  the  above  definition 
the  static  mechanism  is  assumed  to  hold  if  the  peak  overpressure 
Pm,  as  shown  for  a typical  blast  loading  of  Figure  6,  falls  in 
the  range,  Ps  ^ Pm  ^ 3pg.  If  pjj^  > 3pg  the  traveling  plastic 
hinge  mechanism  is  assumed  to  operate. 


Pressure 

P 


Time  t 


Figure  6.  Typical  Blast  Pressure  Time  Curve 
for  a Fuel-Air-Explosion 


_ For  a fixed  end  beam  where  end  rotation  is  by  a plastic 
additional  plastic  hinge  also  forms  at  the  center 
and  the  deflected  shape  is  the  same  as  Figure  4.  The  veloc- 
ity of  the  midpoint  is  found  to  be 


V = 3(T-pgt)/2m 


where  m is  the  mass  per  unit  length  and  T is  the  impulse  per 
unit  length  at  the  time  t.  The  time  t2  where  motion  ceases 


may  be  determined  by  setting  Equation  (7)  to  zero  and  integrat- 
ing the  resulting  equation, 


r ^2 

J p(t)dt  = Pst2  . 


where  I2  is  the  value  of  the  impulse  at  t2  and  for  some  cases 
may  be  determined  by  integration.  The  value  of  t2  may  be 
interpreted  graphically,  as  shown  on  page  32  of  Reference  3, 
for  those  cases  where  an  analytical  expression  of  the  blast 
loading  is  not  known.  Knowing  the  time  t2,  the  final  midpoint 
deflection  may  be  determined  from  the  integration  of  the  veloc- 
ity expression  over  the  time  limits  zero  to  t2*  This  integra- 
tion yields  the  maximum  center  point  deflection 


y(L,t2)  = (/  ^Idt  - Pst^/2] 


The  equations  governing  the  static  mechanism  for  a simply 
supported  beam  are  the  same  as  those  of  the  fixed  end  beam 
except  the  static  collapse  load  for  this  case  is  given  by 
Equation  (6). 

Equations  of  motion  for  a fixed  end  beam  with  the  travel- 
ing plastic  hinge  mechanism,  as  shown  in  Figure  7 and  determined 
by  Abrahams on  et  al^3)^  given  as 


mv  = p , 


Xh  < X £ L 


(10) 


<j^^(I)/3  = pXjj^/2  - 


0 < X < Xt 


(11) 


where  p is  the  general  expression  for  pressure  as  a function  of 
time,  xj^  is  the  position  of  the  plastic  hinge,  and  w is  the 
angular  velocity  of  the  segment  AB  of  Figure  7b.  Using  a con- 
tinuity condition  at  the  hinge  point  B of  Figure  7a  the  velocity 
of  the  underformed  portion  BC  may  be  expressed  as 


V = wxh 


(12) 


With  the  use  of  Equations  (11)^ and  (12)  the  position  of  the  plas- 
tic hinge  xj^  and  its  velocity  Xj^  are  determined,  respectively,  as 


Xh^  = i2Mpt/I 


(13) 
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Figure  7.  Traveling  Plastic  Hinge  Mechanism 
and  Free  Body  Diagram 


Xh  = 6Mpcr-pt)/r2xj^ 


(14) 


The  time  ti  when  the  plastic  hinges  arrive  at  the  center  of 
the  beam  is  determined  by  center  of 


^1  = \(t)dt  = 3p3ti  , 


(15) 


= ri/m 

^ (16) 


'1  = -ti 


'<‘'.ti)  = i /\dt  . 

m Jn 


(17) 


Motion  , of  the  beam  beyond  this  time  +■ 

static  mechanism  and  again  U th^ve^  't 
set  to  zero,  the  time  t ^ expression  (7)  is 

tion  (8).  ihe  finarop^I  ^^en  motion  ceases,  is  given  by  Equa- 
the  expression  Point  deflection  is  then  given  by 


y(L.t,)  = t ± [r^idt-^  (t^2.,2)J 


(18) 


■ *^‘.‘<B^*RKrX'iirrav 


■ .Lr  ■ ■ 1.-:^'’-^  ' 


For  a simply  supported  beam  the  equations  of  motion  are  the 
same  as  those  of  Equations  (10)  and  (11)  except  that  -2Mp  is 
replaced  by  -Mp  due  to  the  operation  of  only  one  plastic  hinge 
at  the  center  of  the  beam. 


As  shown  in  Figure  7 the  preceding  beam  analysis  is  based 
on  a beam  free  of  end  constraints.  This  condition  was  taken 
into  account  in  the  design  of  an  experiment  described  in  Refer- 
ence 3 to  evaluate  the  preceding  analysis.  The  experimental 
arrangement  was  such  that  the  ends  of  the  beam  were  fixed  to 
sliding  supports  which  permitted  the  beam  ends  to  move  toward 
each  other  and  prevented  a buildup  of  tension  load.  Uniformly 
distributed  impulse  values  of  0.06  to  0.146  Ib-sec/in  were 
applied  to  2024-T4  aluminum  beams  by  use  of  thin  strips  of  sheet 
explosive.  The  ratios  of  experimental  centerpoint  deflection 
to  the  theoretical  values  ranged  between  0.402  and  0.652  for 
the  pinned  ends  beams  and  0.517  and  0.752  for  the  fixed  ends 
beams  with  the  higher  ratios  corresponding  to  the  higher  impulse 
loads.  Although  the  analysis  over-predicts  the  centerpoint, 
and  may  not  be  strictly  applicable  to  beam  analysis  where  ten- 
sion is  a factor,  it  does  serve  as  an  important  guide  to  analy- 
tical methods  using  the  plastic  hinge  mechanism. 


Beam  response  to  impulsive  loads  neglecting  axial  restraint 
was  also  studied  by  Symonds'^O^  and  in  a later  paper^H'^  he  in- 
cluded the  effect  of  axial  restraint.  In  this  study  (Reference 
11)  the  assumption  of  traveling  plastic  hinges  was  also  used. 
However,  in  this  case  the  central  portion  of  the  beam  was  given 
an  initial  velocity  Vg.  The  assumption  of  an  initial  velocity 
simplified  the  analysis  and  the  final  shape  of  the  beam  may  be 
found  somewhat  easier  than  the  preceding  analysis.  For  the 
fixed  end  beam  without  axial  restraint  the  final  center  point 
deflection  y(L,t2)  and  rotation  0(L,to)  along  with  support 
rotation  0(O,t2)  are  found  to  be 


y(L,t2)  = 


3(2Mp) 


(19) 


0(L,t2)  = 


mLV^2 


6(2Mp) 


(20) 


0(0, t2)  = 


2(2Mp) 


(21) 


■a 


simply  supported  case  without  axial  constraint  the 
(19)!  (20raS^(21)^  replaced  by  Mp  in  Equations 

;,n  constraint  was  applied  by  SymondsfU)  by  adding 

fSrce  definefas  ''''  M.  A plastic  axifl 


N. 


■dA 


(22) 


IS  assumed  to  be  in  effect  when  the  entire  beam  rross  section 


JL  + Ji  = 1 

Mp  Np2 


(23) 


oocuJT?l^e?er  bf ised^in  IpplJing'^EqSa^iS™ 23  ) 


N 


£ = oli 


Mp  tjj 


N, 


(24) 


terms  of  strain  rate  e and  curvature  lO  was  also 
Fi|S?e  8°fo^  I fixL’'and°S®®  analysis.  Using  the  diagrL  of 
acceleration  of  se^en^^B^^LSrt^JrpoJ^t^fir 

|mx^w  = -2M  - Ny  . 


(25) 


Figure  8.  Beam  Response  Mechanism  Including  Axial  Constraint 
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It  can  be  shown  from  Equation  (24)  that  the  relation  between 
linear  strain  and  curvature  becomes 


N _ Np 

Np  " 2Mp 


(26) 


Using  a velocity  continuity  such  as  Equation  (12)  and  Equations 
(23),  (25),  and  (26),  the  values  of  M,  N,  w,  and  Xj^  may  be  deter- 
mined. Simultaneous  solution  of  these  equations  yields  a third 
degree  equation  in  terms  of  the  final  center  point  deflection 
which  can  be  given  as 


1 ^ 

"NpL 

2 

"yf* 

L 12 

_(2Mp)_ 

_L_ 

^ _ 1 mLVo^ 
3 (2Mp) 


(27) 


where  yf  is  the  final  center  point  deflection  defined  previously 
as  yf  = y(L,t2).  However,  the  final  deformed  shape,  strain,  and 
curvature  functions  require  a numerical  integration  procedure. 
Again,  the  equation  for  final  deflection  of  a simply  supported 
beam  is  obtained  by  replacing  2Mp  with  Mp  of  Equation  (27).  The 
foregoing  analysis  is  based  on  small  deflections  and  is  limited 
to  deflections  on  the  order  of  the  depth  of  the  beam;  hovzever, 
the  important  point  found  here  is  that  the  deflection  based  on 
a model  with  no  axial  constraint  may  be  as  high  as  four  times 
the  deflection  based  on  a model  including  axial  constraint. 
Symonds^ll^  also  discusses  a transition  of  the  beam  response 
of  bending  with  axial  restraint  to  that  of  a plastic  string  with 
end  restraint.  The  assumption  allows  for  calculation  of  deflec- 
tions on  the  order  of  1 to  10  times  the  beam  depth  and  these 
deflections  agree  quite  well  with  results  presented  in  a later 
paper  by  Symonds  and  Jones which  includes  corrections  for 
strain  rate  sensitivity.  However,  the  deflections  in  these 
references  are  still  quite  small  in  comparison  to  gross  deflec- 
tions which  would  be  expected  to  cause  failure. 

The  majority  of  the  previous  researchers  have  avoided  the 
treatment  of  beams  for  loads  intense  enough  to  cause  failure 
from  a structural  standpoint.  This  is  not  to  say  that  intense 
impulsive  loadings  have  not  been  studied  but  much  of  the  work 
in  this  area  has  been  concerned  with,  basic  material  problems 
rather  than  overall  structural  response.  However,  tearing  and 
shear  failures  in  explosively  loaded  clamped  beams  were  studied 
by  Menkes  and  Opat^l^^  and  a failure  criterion  based  critical 
particle  velocity  is  given  by  Sewell  and  Kinney (1^^. 
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Slo^JIrsf caJiLTirIng!  ^As^^hfi^^LsS'^y^cl 

increases,  overlapping  of  modes  ? e„a  a J!  • * impulse 

, Upon  examination  of  results  of  this  wnr^v(13)  ^ 
points  were  observed.  work  several  mam 

Tom  shape  of  the  beams  subjected  to  the 

low  impulse  intensities  indicates  a static  ollstlo 

d^SmalSrprocesL"^"  final 

^ ^^he  result  of  tensile 

thrsamJ  Jar?,M®  2 threshold  is  dependent 
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3 3 (shear  failure)  is  dependent  on  the  initial 
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tensile  critical  particle  velocitv(5^y  ’ exceeded  the 

new  failure  criterion  ^al  Ilp^ied  to  plSL  and’”r®'^^%^- 
When  the  positive  oressu-rp  p 4-u  beam  elements 

than  one-quarter  of  the  cha^po+-pi-)°  +-*^^  blast  wave  was  less 
element.  However  it  wpc;  ^ fistic  period  of  the  structural 

the  Shear  cSlIcai  pL??Lf  ve?olS^"L*’X  1^"  tha? 

16,  gave  better  predictions  than  +ha+  I defined  in  Reference 
particle  velocitrin  exneSLe^?.  “at  of  the  tensile  critical 
panels.  A comparison  of  the  shJaJ'’fa^T^"®  failure  of  aluminum 
reported  in  Reference  13  Sd  ?hS?  of?  impulse 

critical  particle  velocity  is  given  irF!g:re“""®TjfsSir 
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Figure  9.  Impulse  Required  for  Shear  Failure  Mode  of  Beams 


critical  impulse  curve  of  Figure  9 is  based  on  a shear  critical 
particle  velocity 16 ) given  as  ^ 

2 

(29) 

where  is  the  tensile  critical  particle  velocity,  is 

the  ultimate  tensile  stress  and  p is  material  density.  The 
corresponding  critical  impulse  for  shear  failure  is  then 

~ ^^CR^  mh  , (30) 

s s 

where  h is  the  beam  depth.  The  tensile  critical  particle 
velocity  for  some  aluminum  alloys  is  approximately  220  ft/sec. 
cy  Rinehart . 


(vcr)^  = (vcr)^ 


'u 


<Vcr>, 


- 1 


Based  on  the  precedii^  discussion  of  previous  work  it  is 
apparent  a need  exists  for  some  analysis  to  predict  gross  de~ 
flections  for  mild  impulse  intensities  found  in  fuel-air  explo- 
sions . The  next  subsection  is  devoted  to  a strength  of  materials 
approach  for  predicting  beam  response  in  terms  of  plastic  hinges 
and  axial  restraint  with  la??ge  deflections. 
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1.4  Response  of  Beams  Subjected  to  a Fuel^Air  Explosion 

For  loading  functions  of  fuel-air  exnloQione  ‘j.* 

pressure  phase  is  rather  long  (milusLoS  ringed  Sd 
Of  the  pressure  time  function  becomes  important  It  irthS  ^ 

air  explosions  analyzing  beam  response  to  fuel- 

».ent  of  a bean.  Vesponsf mod^^whe^subjeStef ?o"rJef? 

ri^-L  d!^?; 

arrf?xed"agalnfr2rmSvirn?  SwLd”eacrot?er^"  ^otr^he®®"" 

fxxed  end  beam,  allowing  no  rotation  except 

Snges?  »>'  natHjSf  ’ 

energy  is  much  larger  thaA  the  elastic  eneiJ^yf 

plastifhSgeffo™  afthf«n?erand‘*JnS’‘l?  ?hfb°a 

nleSed^fS/rfagrlngian^forLlaholl^rttfequa^ 

may  be  written  Il<;o  t • equation  of  motion 

y wriTTen.  Use  of  the  Lagrangian  method  is  justified  by 
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Figure  10.  Free  Body  Diagram  for  a Static 
Analysis  of  a Beam  Fixed  at  Both  Ends 

coLdinatf  system^and^by^writillg^th^  undeformed 

terms  of  the  plastic  enLgy  ® dissipative  functions  in 


The  following  symbols  will  be  used  in  the  derivation  of 
the  equation  of  motion  for  the  half  span  of  Figure  6 . 

A cross  section  area  of  beam 


L half  span  length 

Mp  plastic  moment 

My  yield  moment 

moment  of  inertia  of  half  span  about  the  end, 

c .. 

■ 3 

I(t)  impulse-time  function  (psi-msec) 

p(t)  pressure-time  loading  function  (psi) 

V Velocity  of  midspan  (in/sec) 

w width  of  beam  (in) 

6 deflection  of  midspan  (in) 

A deflection  along  beam  longitudinal  axis  (in) 

ay  yield  stress  (psi) 

p mass  density 

0 angular  rotation  of  beam  (rad) 

x,y  coordinate  axes  along  original  beam  position. 

The  kinetic  energy  KE  may  be  written  assuming  a uniform 
linear  velocity  variation  in  the  direction  of  longitudinal 
axis  of  the  half  span  and  is  expressed  as 


= . I r if)  . 


(31) 


Making  use  of  the  moment  of  inertia  given  above,  Equation  (31) 
reduces  to 


KE  = (e^L^  + A^)  . 

6 


(32) 
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and 


Now  A and  A may  be  expressed  as 

A = (sec0-l)L 
• • 

A = (Ltan0sec0)0 

Using  Equations  (32)  and  34)  the  kinetic  energy  becomes 

3 


(33) 

(34) 


KE  = (l+tan2©sec20)02  . 

6 


(35) 


The  plastic  energy  PE  is  now  written  as 

PE  = 2Mp0  + AcyA  , (30) 

and  using  Equation  (33)  the  plastic  energy  may  be  written  as 

PE  = 2Mp0  + (sec0-l) )0yAL  . (37) 

pe  generalized  force  determined  from  the  applied  pressure 
force  IS  expressed  as 

(5)  p(t)wdx 


Q0 


iU 


(38) 


Assuming  no  shortening  of  the  span  then 
6 = Ltan0  , 

and  the  generalized  force  becomes 


(39) 


T = sec2e  . 


(40) 


Defining  the  Lagrangian  L and  the  Lagrange  equation  as 
L = KE  - PE 

and 


(41) 


(42) 


1 


i 


j i 
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(l+tan^Gsec^©)© 

3 


+ (sec^0tan0+tan3©sec^©)©^ 

O 


I 


I 


+ 


2Mp  + ayAL( tan0sec0) 


_ p(t)wL^^  2r> 
= sec^Q 

2 


(43) 


Equation  (43)  is  valid  only  for  positive  values  of  ©.  The 
basic  assumption  is  that  a plastic  hinge  is  in  effect  at  the 
beginning  of  the  beam  motion.  This  means  that  at  t = 0 the 
moment  resulting  from  the  initial  blast  loading  must  be  large 
enough  to  overcome  the  internal  resisting  force  2Mp.  For  Equa- 
tion (43)  a necessary  condition  to  insyre  positive  initial 
values  of  0 is  that  initial  values  of  0 be  positive.  This 
criterion  may  be  met  by  the  requirement  that 


Pmax  1 4Mp/wL2 

or  (44) 

Pmax  1 4f0yz/wL2  . 

The  assumption  is  also  made  that  the  maximum  value  of  0 is 
reached  in  the  interval  0 < t < t.  If  this  assumption  is  not 
true  then  an  equation  with  the  right  hand  side  or  forcing 
function  equal  to  zero  would  hold  for  this  case.  This  case 
is  omitted  in  that  either  maximum  deflection  or  failure  is 
expected  to  occur  within  the  interval  0 < t < t . 

The  strain  induced  in  the  beam  due  to  this  deflection  is 
composed  of  an  average  axial  strain  obtained  from  Equation  (33) 
and  defined  as 

^a  = ^ = (sec0-l)  , (45) 

and  a strain  £0  at  the  outermost  fiber  at  the  plastic  hinge 
points . An  estimate  of  the  magnitude  of  £3  will  be  made  from 
the  following  simple  analysis. 


I 
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In  the  preceding  analysis  a plastic  hinge  was  assumed 
to  be  concentrated  at  a point;  however,  in  reality  the  hinge 
point  is  the  position  where  the  yield  stress  first  extends 
over  the  entire  depth  of  the  cross  section.  Extending  out 
from  either  side  of  the  hinge  point  is  a plastic  zone  of  total 
length  Lp  as  shown  in  Figure  11a.  This  zone  is  essentially  a 
volume  of  the  beam  in  which  the  internal  stress  is  equal  to 
the  yield  stress. 
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Plastic  Hinge 
Location 


Figure  11.  Schematic  of  Plastic  Zone 
Reference'^9  8“®" 


Lp  = 2L(l-l/f)  , 

beam  span  and  f is  the  beam  shape  factor, 
fa^*’  length  corresponding  to  the  angle  20 

-0  = I . 


(46) 


(47) 


The  x’adius  of  curvature  R is  then 


R = 


(48) 


Equations  (46)  and  (48)  into  Equation  (47) 
gj ves  the  strain  due  to  rotation  as 


£0  = 


L(l-l/f) 


(49) 


Defining  strain  e as 


the  total  strain  at  the  outermost  fiber  of  the  midpoint  and 
fixed  end  of  the  beam  may  be  expressed  in  terms  of  0 as 


e(0)  = e(L)  = (sec0-l)  + 


C0 


L(l-l/f) 


(51) 


The  basic  assumption  in  this  case  is  that  the  fixed  end  beam 
is  fixed  against  rotation  and  translation,  whereas  the  simple 
supported  beam  is  fixed  against  translation  only. 


The  equation  of  motion  for  the  simple  supported  beam  will 
be  the  same  as  Equation  (43)  except  the  third  term  of  the  right- 
hand  side  will  be  replaced  by  Mp.  The  resulting  strain  for  the 
midpoint  of  a simple  supported  beam  will  be  the  same  as  Equation 
(51);  however,  for  the  support  point  of  a simple  supported  beam 
the  strain  will  be  only  (sec0-l),  due  to  the  operation  of  a 
natural  hinge  at  the  support. 


Solution  of  Equation  (43)  was  obtained  using  a fourth  order 
Runge-Kutta  analog  simulation  program.  Solutions  were  determined 
for  6061-T6  aluminum  beams  of  the  same  size  as  those  of  Refer- 
ence 13.  Maximum  strain  values  were  determined  using  Equation 
(51)  and  any  strain  exceeding  0.15  was  assumed  sufficient  to 
cause  failure.  Results  of  these  calculations  along  with  the 
experimental  and  theoretical  results  of  Reference  13  are  listed 
in  Table  2.  Experimentally  determined  pressure  time  functions 
of  Reference  6 were  used  as  an  applied  load  p(t)  of  Equation 
(43) . 


Applying  the  same  general  approach  to  that  of  a traveling 
plastic  hinge  as  shown  in  Figure  12,  the  kinetic  energy  may  now 


cyA  Mp 


(a) 


Figure  12.  Schematic  of  Traveling  Plastic  Hinge  Mechanism 
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TABLE  2.  BEAM  RESPONSE  TO  IMPULSIVE  LOADS 


Impulse 

Ktaps 


Length  Thickness 
(Inches)  (Inches) 


^sntral  Deflection 

Reterence  13 

Exp  Equat 

(Inches)  (Inches)  (Inc 


10.9 

8 

0.187 

0.99 

0,94 

15.3a 

8 

0.187 

+ 

+ 

17.8 

8 

0.187 

1.56 

1.44 

10.9 

8 

0.250 

0.54 

0.50 

15.3a 

8 

0.250 

+ 

+ 

17.8 

8 

0.250 

1.07 

1.50 

15.3a 

8 

0.375 

+ 

+ 

17.8 

8 

0.375 

0.72 

0.50 

10.9 

4 

0.187 

0.44 

0.25 

15.3a 

4 

0.187 

+ 

+ 

17.8 

4 

0.187 

0.77 

0.69 

10.9 

4 

0.250 

0.32 

0.31 

15.3^ 

4 

0.250 

+ 

+ 

17.8 

4 

0.250 

0.53 

0.44 

15.3a 

4 

0.375 

+ 

+ 

17.8 

4 

0.375 

0.48 

0.18 

Equation  43 
(Inches ) 


+ 

0.769 

+ 

+ 

0.475 

+ 

0.154 

+ 

+ 

* 

+ 

+ 

ft 

+ 

ft 

+ 


*No  response 

+No  calculation  or  value 


p(t)  = 830(l-t/.0015)  -t/.0015 


be  expressed  in  terms  of  the  angle  0 and  the  plastic  hinge  posi- 
tion x^.  Assuming  that  the  central  portion  of  the  beam  BC  of 
Figure  12a  remains  rigid  and  resists  the  loads  oyA  and  Mp  with- 
out deformation,  the  kinetic  energy  KE  and  the  plastic  energy  PE 
may  be  defined  in  a manner  similar  to  Equations  (31),  (32),  (33), 
(34),  (35),  (36)  and  (37).  The  plastic  energy  in  this  case  is 
the  same  as  Equation  (37)  except  L is  replaced  with  xv.  An  ex- 
pression for  the  kinetic  energy  must  now  include  a term  for  the 
central  portion  of  the  beam  as  well  as  the  portion  of  length 
Xh*  The  expressions  for  the  energy  terms  are  found  to  be 

KE  = [xj^d+tan^esec^G)  + 3(L-Xh)sec‘^0]0^ 

+ pAxj^  [ (L-xj^)sec‘^0tan0]Xj^0  (52) 

+ qA  [ (L-Xj^)tan^0]x^ 

and  PE  = 2Mp0  + (secG-DayAXj^  . (53) 

Using  the  continuity  of  deflection  at  point  B of  Figure  12, 

6 = xj^tanG,  the  force  component  expression  for  the  two  coordi- 
nates and  X],  may  be  determined.  Using  the  total  work  done  by 
the  pressure  force  the  force  components  are 

Q0  = xj^  (^“^)  wp(t)sec20 

and  (54) 

Qxh  - (L-Xj^)wp(t )tan0  . • 

Using  Equations  (41),  (42),  (52),  (53)  and  (54)  the  equations 
of  motion  may  be  determined  and  are  given  as 

pA  [(L-Xh)tan20]xh  + pAxj^  [ (L-xj^)sec20tan0]0 

+ 2pA  [ (L-X2^)sec^0tan0]0xj^ 

+ pAxj^  j^(2L-|x2^)  sec^Gtan^G  + ^^ec'^G-^J©^  (55) 

- ^ [tan^0]x^  + aYA(secG-l) 


j 

I 


I 

i 


= p(t  )w(L-Xi,^ ) tanG 


[Y(l+tan20sec20)  + (L-Xh)sec‘*0je 
+ pAxj^  C(L-Xh)sec20tan0]xjj 

+ pAx'n  Cxh(l+tan20sec20)  + (2L-3xj^)sec‘*0]0Xh  (56) 

+ pAxg  C2(L+Xj^)sec‘*0tar9  + ^(tan0sec‘^0+tan30sec20)]e2 

” pAxh  [sec20tan0]x2  + 2Mp  + OYAxhSec0tan0 

= p(t)wxh  (l~)  * 

No  attempt  was  made  to  obtain  a solution  to  these  equations. 
1«5  Conclusions 
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material  constitutive  relation  may  be  used  when  the  nlastic 

?.r?» 

sis  may  be  in  order.  Tests  of  beams  and  plates  stiffened  with^ 
beam  elements  subjected  to  actual  fuel-air  exploLJns  (bL 
periments)  are  recommended.  ^ ^bag  ex- 


SECTION  II 

FUEL-AIR  EXPLOSION  (FAE)  IN  CONFINED  STRUCTURES 


2 . 1 Introduction 

Before  considering  the  effects  of  fuel-air  explosions  in 
large  closed  or  confined  areas  a discussion  of  the  fuel-air 
explosion  in  unconfined  areas  is  in  order.  For  a given  fuel-to- 
air  mixture  many  of  the  characteristics  found  in  a fuel-air  ex- 
plosion are  applicable  whether  the  gas  cloud  is  in  a large 
container  or  open  to  the  atmosphere. 

2.2  Fuel-Air-Explosion 

A fuel-air  explosion  is  a rather  complicated  phenomena; 
however,  in  broad  simple  terms  of  explanation:  given  a proper 

fuel-to-air  mixture  and  a small  high  energy  secondary  explosion 
at  some  point  or  line  in  the  mixture,  a detonation  wave  is 
established  and  moves  through  the  fuel-air  mixture  at  some  con- 
stant velocity.  Under  the  proper  conditions  the  wave  behaves 
like  a strong  blast  wave  but  later  it  may  be  described  as  a 
Chapman-Jouguet^^^^  detonation  wave.  The  velocity,  pressure, 
temperature,  and  density  changes  across  the  wave  are  all  step 
functions  as  in  the  case  of  a strong  shock  wave  in  air;  however, 
in  the  case  of  the  detonation  wave  they  are  all  constant  with 
time  and  depend  only  on  the  properties  of  the  fuel-air  medium 
ahead  of  the  shock.  The  basic  difference  between  the  detona- 
tion wave  and  the  regular  shock  is  that  the  energy  equation 
for  the  detonation  wave  has  a term  which  represents  the  heat 
liberated  in  the  reaction  zone  and  is  added  to  the  energy  flux 
passing  through  the  reaction  zone (18).  This  heat  addition  com- 
ing at  the  wave  front  causes  a larger  temperature  jump  than 
that  of  a regular  shock.  In  turn,  this  increased  temperature 
increases  the  acoustic  velocity  and  decreases  the  particle 
velocity  behind  the  front. 

The  wave  may  be  planar,  cylindrical  or  spherical,  depending 
on  the  type  of  secondary  detonation.  The  static  or  side-on 
pressure  jump  across  the  wave  is  the  same  in  all  three  and  is 
constant  with  time.  Only  the  shape  of  the  wave  behind  the  front 
is  changed  with  time  or  distance  from  the  secondary  explosion 
source  and  is  highly  dependent  on  distance  to  the  nearest  re- 
lief surface.  For  a strong  detonation  wave  the  static  or  side- 
on  pressure  ratio  of  pressure  across  the  wave  front  is  given 
by 
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(57) 


Pi  . pqD^  . i-»-yoD^ 
PO  YO‘^1  Yi+1 


where 

Pq  - static  pressure  in  medium  ahead  of  the  wave  (psi) 

P]^  - static  pressure  at  wave  front  (psi) 

PO  - density  in  medium  ahead  of  the  wave  (Ib-sec^/in**) 

D - wave  front  velocity  (in/sec) 

Yo  - ratio  of  specific  heats  in  front  of  detonation  wave 

Yi  - ratio  of  specific  heats  behind  detonation  wave. 

The  reflected  pressure  Pr  exerted  on  a rigid  wall  given  in 
Reference  18  for  a condensed  explosive  and  verified  experimen- 
tally by  Ross  and  Strickland ^ 19'  for  a fuel-air  mixture  may  be 
expressed  in  terms  of  the  side-on  pressure  as 

Pt^Syi^i  + ✓i7yi^2Yi-H  ^ 

Pi  4yi 

As  shown  in  Equation  (58)  the  ratio  of  reflected  pressure  to 
the  side-on  pressure  is  dependent  only  on  the  ratio  of  specific 
heats  and  is  fairly  insensitive  to  changes  in  values  of  y above 
.unity.  The  value  of  pp/p]^  is  2.6,  2.4  and  2.3  for  y values  of 
1,  3 and  »,  respectively.  This  indicates  that  the  reflected 
pressure  will  be  about  2.4  times  the  static  pressure  for  almost 
any  fuel-air  mixture.  Further  observations  for  fuel-air  explo- 
sions^!'' show  that  for  most  fuel-air  mixtures  the  detonation 
velocities  are  approximately  5000  to  6000  FPS.  This  means  that 
for  almost  any  fuel-air  explosion  the  reflected  pressure  is 
nearly  a constant  and  will  be  much  lower  than  that  of  a reflect- 
ed shock  wave  in  aii;  with  the  same  static  pressure  rise. 

The  static  pressure  rise  and  corresponding  reflected  pres- 
sure siyen  by  Equatiohs  (57)  and  (58)  were  verified  experimen- 
tally^ ' for  a MAPP-Air  system.  Impulse  calculations  based 
on  the  experimental  pressure  time  histories  also  check  with 
analytical  predictions  of  Sichel  and  Hu^^O),  Analytical  pre^ 
dictions  of  fuel-air  explosions  are  given  by  Nicholls,  et 
and  the  reader  is  referred  to  these  documents  for  further  dis- 
cussion of  unconfined  fuel-air  explosions  in  the  atmosphere. 
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2.3  Confined  Fuel-Air  Explosions 

In  the  previous  section  a fuel-air  explosion  was  defined 
in  terms  of  a detonation  wave  travsling  through  a fuel-air 
medium.  One  of  the  major  assiunptions  was  that  the  fuel-air 
cloud  was  of  sufficient  size  and  mixture  to  produce  a strong 
detonation  shock.  For  the  confined  case  the  assumption  will 
be  that  the  limits  of  the  confinement  are  large  enough  such 
that  the  same  conditions  are  satisfied  to  produce  a detonation 
wave,  and  the  mechanism  will  be  the  same  as  the  unconfined  case 

In  the  case  of  unconfined  fuel-air  explosions  the  ground 
impulse  and  the  dynamic  oi'  reflected  impulse  are  often  the 
parameters  used  to  describe  the  dairage  capability  of  the  ex- 
plosion. It  is  assumed  that  this  v ill  also  be  true  for  the 
confined  fuel-air  explosion.  In  tie  analysis  of  an  idealized 
unconfined  fuel-air  explosion^22)  it  was  assvimed  that  no  side 
or  top  relief  occurred.  This  is  essentially  the  same  assump- 
tions that  would  be  used  in  the  analysis  of  a confined  fuel- 
air  explosion.  This  analysis  gives  a typical  side-on  pressure 
time  history  or  shape  as  chown  in  Figure  13.  The  effect  of 
confinement  on  a FAE  is  to  maintain  some  high  pressure  for 
some  finite  time  as  evidenced  by  the  plateau  of  the  pressure 
time  curye  shown  in  Figure  13.  In  the  unconfined  case  a relief 
wave  moving  at  the  acoustic  speed  in  the  burned  gas  quickly 
reduces  the  pressure  behind  the  detonation  wave  and  the  result- 
ing pressure  time  curve  is  similiar  to  that  of  Figui’e  6.  For 
a confined  FAE,  the  reflected  pressure  profile  will  be  similar 
in  shape  to  that  of  Figure  13  except  the  peak  overpressure  will 
be  approximately  2.4  times  that  of  p^^.  The  ground  inipulse  and 
dynamic  or  reflected  impulse  as  calculated  in  References  20  to 
22  will  be  applicable  to  fuel-air  explosions  in  buildings,  air- 
craft shelters,  hangers,  etc.  In  such  cases  a detonation  wave 
whose  shape  depends  on  location  and  size  of  the  secondary  ex- 
plosion will  move  through  the  medium  at  a constant  velocity 
and  be  reflected  from  the  walls  and/or  top. 

The  wave  type,  i.e.,  planar,  cylindrical,  spherical,  will 
depend  on  the  general  shape  of  the  building  and  type  and  loca- 
tion of  secondary  explosion.  For  rather  long  buildings  there 
will  be  a tendency  to  form  a planar  wave  as  the  detonation 
front  moves  down  the  length  of  the  structure.  In  this  case 
the  major  damage  would  be  caused  at  one  end  due  to  the  reflec- 
tion of  the  wave.  In  most  buildings  roof  lift  may  be  the  most 
vulnerable  mechanisms  and  in  such  cases  an  initiator  or  second 
event  charge  on  the  floor  would  give  wave  reflection  at  the 
roof  or  overhead.  Due  to  the  rather  small  speak  overpressure 
of  FAE  devices,  it  is  essential  to  try  and  direct  the  wave  in 


such  a direction  as  to  cause  the  most  damage.  Of  course,  this 
may  be  impractical  to  do  and  in  turn  calculation  of  the  exact 
pressure  or  impulse  at  a point  inside  the  structure  may  not  be 

possible  even  if  the  exact  location  of  the  secondary  blast  was 
known. 


Figure  13.  Pressure  of  a Fuel-Air  Explosion 
Versus  Distance  from  the  Secondary  Explosion 


For  point  secondary  explosions  or  initiators  a spherical 
wave  will  be  established  and  if  the  geometry  of  the  reflect- 
ing surface  is  given  a reasonable  estimate  of  the  impulse  may 
be  determined  based  on  calculations  of  Reference  20. 


2.4  Conclusions 

K®fl6cted  pressures  of  fuel— air  explosions  from  rigid 
walls  will  reach  a maximum  of  approximately  2.4  times  the 
static  pressure  rise  of  the  detonation  wave.  The  same  detona- 
tion wave  mechanism  of  unconfined  fuel-air  explosions  without 
edge  relief  will  be  operative  for  fuel-air  explosions  in  large 
unpressurized  containers.  As  in  the  case  of  unconfined  fuel- 
air  explosions  the  main  damage  mechanism  for  large  structures 
will  be  the  reflection  of  a detonation  from  the  walls,  top,  or 
end  of  the  confining  structure.  For  confined  explosions  where 
occurs  pressure  relief  will  occur  upon  venting  and  the 
effect  of  this  pressure  relief  should  be  taken  into  account 
when  calculating  the  impulse  intensity  or  positive  pressure 
phase  of  the  detonation  wave. 


SECTION  III 


STRUCTURAL  DEGRADATION  OF  CONCRETE  BY  IMPULSIVE  LOADS 


3 . 1 Introduction 
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Section  3.2  considers  the  modeling  of  surface  and  near 
surface  shots  in  rock  and  concrete.  Most  of  the  information 
given  is  for  rock  rather  than  concrete,  because  that  is  what 
is  available.  It  is  believed  that  the  same  procedures  can 
be  applied  to  concrete  if  sufficient  experimental  data  can 
be  obtained. 


3 . 2 Empirical  Cratering  Formulas  for  Surface  and  Near  Surfa-^e 
Shots 

3.21  Preliminary  Comments 

This  part  of  the  investigation  addressed  the  question  of 
whether  it  is  possible  to  predict  the  crater  volumes  and  dimen- 
sions produced  by  surface  explosions  on  massive  concrete  blocks 
or  slabs  by  the  use  of  simple  empirical  formulas.  Not  enough 
experimental  data  on  concrete  was  found  to  test  the  empirical 
formulas,  but  there  are  several  reports  on  cratering  in  rock 
by  HE,  both  from  blasting  research  related  to  mining  and  ex- 
cavation and  from  extensive  studies  of  the  modeling  of  nuclear 
explosive  loading  by  HE  experiments.  These  results  lead  to 
the  conclusion  that  cratering  in  concrete  could  be  modeled  in 
the  same  manner. 

The  empirical  formulas  have  quite  successfully  represented 
true  crater  dimensions  as  a function  of  charge  weight  over  a 
range  of  one  pound  to  100  tons,  provided  that  the  material  was 
reasonably  homogeneous  (no  large  faults  or  layers  of  different 
material)  and  that  the  surface  contact  conditions  and  charge 
shape  and  detonation  initiation  position  were  maintained  the 
same.  The  nuclear  explosive  modeling  encountered  difficulties 
with  extension  up  to  the  megaton  range (33),  but  the  model 
studies  themselves  showed  remarkable  consistency. 

It  is  recognized  that  successful  empirical  representation 
of  controlled-contact-condition  surface-burst  loadings  will  be 
of  limited  value  for  practical  prediction  of  damage  from  bombs. 
Surface  explosive  attack  is  not  an  efficient  way  to  remove 
material.^  A buried  charge  is  much  more  effective,  and  there 
is  an  optimum  depth  of  burial  for  a given  charge  and  a given 
material.  To  maximize  material  removal  by  cratering  one  would 
seek  to  have  the  bomb  penetrate  to  the  optimum  depth  before 
exploding.  This  introduces  a host  of  problems  concerned  with 
projectile  design  and  fuzing  and  shaped  charge  design  for  effec- 
tive penetration  which  are  not  considered  in  this  report. 


It? 


I Waterways  Experxment  Station  compendium( 33 ) com- 
appear  sufficient  to  predict  bomb/shell 

formed  by  bare  charges  equivalent 
yield  to  the  bomb  or  shell  filler  material,  even  though  it 

1°«  in  rupturing  thfcase 

sive  fillSr  of  the  fragments.  The  explo- 

sive  filler  IS  usually  cast  in  a cylindrical  configuration,  and 

head^^tbpS^.^^  usually  initiated  in  the  nose  or  tail  of  the  war- 
head, these  considerations  also  fail  to  answer  completely  the 
questions  raised  over  differences  in  craters.  ^ ^ 

I'nH-Jr,  Study  which  has  been  done  on  this  subject 

deeply  buried  bombs  produce  craters  smaller  in 
radius  but  larger  in  depth  than  do  comparable  bare  charges 

ho2t  be  sufficient  data  on  shallowly  buried 

leJiS  QbPi?®  ^^{perience  with  mortar  and  artil- 

indLpbPc^^^  contain  a proportionately  heavier  casing, 

ib?iS  near  surface  (fuze  quick)  bursts  form  craters 

somewhat  larger  than  those  of  bare  charges.  Unfor- 

■'V 

tonatS:  b«r::irSocii:^ier  it; 

example.  White  (Reference  26,  Chapter  15)  mentions  tests  on 
concrete  plates  using  demolition  blocks  on  end  and  detonated 

The  total  ?Spuisf??ans- 

tion  ft  thP^?nn^^JVJ?  ^ percent  greater  for  detona- 

initiation. ThTtele  lllovT  ^ 
?*  relationship  ? = 1.5J  between  crater  volume  V 

(in  ) and  total  impulse  J (lb-sec)  for  normal  concrete  (com- 
pressive strength  about  4000  psi)  for  charges  on  enrf  ^ 

at  the  end  away  from  the  slab^  but  warS^thlt “hl  fe^ious 
scattering  of  the  data  indicates  that  other  effects  are  import- 

mentioned  was  even  for  fairly  well  controlled 
so  wen?"'-  conditions  woul?  nit  brSoS“®" 


of  further  emphasized  the  importance 

°f  contact  conditions.  For  example,  Haas  and  Rinehart(39)  oeoem 

laboratory  studies  of  the  coupling  from  smiu  IIIIIliLd  anr*^  ” 
partially  confined  cylindrical  charges  (several  grams)  to  uni 


(a)  Length  and  diameter  of  the  charge 

(b)  Obliquely  oriented  charge 
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(c)  Partial  confinement  (i.e.,  partial  insertion  into  a 
hole  drilled  in  the  rock) 

(d)  Attenuation  with  distance  in  the  rock,  and 

(e)  Coupling  through  layers  of  water,  drilling  mud,  rock, 
aluminum,  plexiglass,  vermiculite,  sand,  foam  rubber, 
or  air. 

The  effects  of  diameter  (or  contact  area) , confinement  and  coup- 
ling material  were  especially  noteworthy.  The  cylindrical  charges 
were  formed  by  stacking  discs  of  DuPont  EL-506A. 

The  peak  stress  at  2-inch  depth  for  0.5-inch-diameter  charges 
increased  from  26,500  psi  for  a 0.25-inch-long  charge  to  34,000 
psi  for  1.5-inch-long  charges,  with  no  further  increase  for  3- 
inch  lengths.  This  is  not  a large  increase,  considering  that  the 
1.5-inch-long  charge  weighed  six  times  as  much  as  the  0.25-inch- 
long  charge.  The _ constancy  of  the  peak  stress  for  charge  lengths 
greater  than  1.5-inch  was  explained  by  Haas  and  Rinehart  from 
considerations  of  detonation  head  theory.  For  an  unconfined 
cylindrical  charge  detonated  from  one  end,  a stable  detonation 
head  is  formed  by  the  time  the  detonation  front  has  propagated 
approximately  three  charge  diameters.  Making  the  charge  longer 
(for  a fixed  diameter)  does  not  increase  the  peak  stress 
propagated. 

A measure  of  the  attenuation  of  the  peak  stress  with  distance 
of  travel  was  obtained  by  varying  the  slab  thickness  and  measuring 
the  momentum  per  unit  area  transmitted  to  a small  pellet  at  the 
back  face  of  the  slab.  For  the  0.5  x 1.5— inch  charge  the  peak 
stress  in  the  rock  dropped  from  271,000  psi  at  0.5-inch-depth  to 
about_80,000  psi  at  1 inch,  34,000  psi  at  2 inches  and  15,200  psi 
at  4 inches.  This  decay  is  considerably  faster  than  the  1/r  decay 
curve  that  would  be  predicted  for  a spherical  elastic  wave  from 
a point  source,  especially  between  0.5  inch  and  1 inch. 

The  peak  stress  transmitted  to  a given  distance  of  2 inches 
(from  a 1- inch- long  charge)  was  strongly  dependent  on  the  charge 
diameter,  increasing  from  8010  psi  for  a 0.27-inch  diameter  to 
61,700  for  a 0.75-inch  diameter.  This  was  attributed  to  two 
factors . The  duration  of  the  applied  pressure  is  greater  for 
the  larger  diameter  (1.4  microseconds  for  0.75  inch  compared  to 
0.5  microsecond  for  0.27  inch).  The  longer  pulse  increases  the 
extent  of  crushing  near  the  charge.  The  shock  wave  peak  stress 
emerging  from  the  crushed  region  was  said  by  Haas  and  Rinehart 
to  be  equal  to  the  dynamic  compressive  strength  of  the  rock.  For 
a larger  crushed  region  there  is  a shorter  attenuation  distance 
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I 
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from  the  point  where  the  peak  stress  is  en„:.i 

ing  strength  to  the  measurement  pSint  ft  5 “'’“f’’" 

charge.  The  second  factor  involifd  L ^hf  a 

SftfSc°f!  """=  =Phfrilff%fffp”fL%"frfhf%5ifc“ 


test  the'typMhffif^fffthf^fS^fSfff^f  explicitly 

Of  the  rock.  If  there  is  in  compressive  strength 

dynamic  crushing  strength  Pvr,r.«o*^  a well-defined 

stress  in  TspSIrical  waL  ^^dial 

gate  out  of  tL  crushed  2oAe,^lince  propa- 

stress  would  produce  further  crishiAe  ^ Thf  o'"  higher 

ant  because,  if:  mg*  The  question  is  import- 


simple  scalifg’^lf„f^“aljf‘^  region  can  be  predicted  by 


Sfshff'rfgifrif"taS®  Tf  from  the 

the  =onstifutive%$u"ftLr^Sfe“r\Kf 


will  arrive  st%hf^baok^ffL°ff^I i®''®!  that 

however,  that  the  crashed  fefiof  if  be  noted, 
the  true  crater  estimated  bv  thf>  om  * • general  identical  to 
in  the  following  subsection  sincp  discussed 
the  large  block!  and  cSiArioo!eS!d  includes  also 
the  initial  compressive 


and  tt  fafbff''dfff:ff:d^fL®flf^  the  0.5  x 1.5-inch  charge 
inch  depth  from  33.800  psi\or  dirfft^ff  ^t>ansmitted  to  the  2- 
psi  for  a 1/8-inch  gap  ^21  Gnn  contact  to  about  24,000 

psi  for  a 1-incS  gapf'^’pff^lf?  yV^ch  gap  and’l2.200 

1.5-inch-long  charge  into  a dm* i insertion  of  the 
almost  doublld  the  fefk  ftfef.  “-5  “®b  deep 

end  of  the  charge?  iffsf  T°  ^ “=hes  from  the 

emphasize  the  problem  of  predictiy'^th°'“^ey°'‘  °bservations 
sions  when  the  contact  conditions  Ire  no? 


j 
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The  importance  of  the  contact  conditions  is  also  shown 
in  the  scaling  of  the  surface  bursts  reported  in  Section  3.2-2. 
It  is  necessary  to  distinguish  between  so-called  true  surface 
bursts , where  the  charge  is  half  buried  with  the  mass  center 
of  the  charge  at  ground  level,  requiring  a prior  hemispherical 
excavation  for  emplacement  of  spherical  charge,  and  the  surface- 
tangent  charge  configuration,  where  the  mass  center  of  a spher- 
ical charge  is  at  a height  of  burst  (HOB)  equal  to  one  charge 
radius  r^ . 


3.2-2  Scaling  Laws  and  Empirical  Formulas  for  Cratering 


If  geometric  similarity  is  maintained  by  the  craters  pro- 
duced in  a given  material  by  charges  of  varying  size,  then  each 
linear  dimension  of  the  crater  is  proportional  to  the  cube  root 
of  the  volume.  If  in  addition  the  crater  volume  is  proportional 
to  the  charge  weight  W,  then  the  crater  dimensions  will  be  pro- 
portional to  Wl/3.  This  cube  root  scaling  law,  predicted  by 
dimensional  analysis  on  the  basis  of  cerxain  simplifying  assump- 
tions, has  been  widely  used  (see,  e.g.,  Lampson,  Reference  26, 
page  114) . 


An  extensive  study  of  craters  from  surface  explosions  and 
scaling  laws,  Vortman'^'^^  indicated  departures  from  the 
scaling.  ^Much  more  data  were  available  from  buried  charge  ex- 
plosions in  earth  and  rock  than  from  surface  bursts,  and,  accord- 
ing to  Vortman,  the  available  experimental  surface-burst  data 
(from  a variety  of  chemical  and  nuclear  explosions)  were  occas- 
ionally contradictory.  If  empirical  laws  of  the  form 


n. 


n. 


ra  = CnW 


da  = C2W 


Va  = C3W 


(59) 


are  fitted  to  the  data  for  apparent  crater  radius  r^,  depth  d^, 
and  volume  V^,  the  data  from  nuclear  and  chemical  explosions  are 
consistent  in  indicating  that  n^  is  larger  than  one-third. 
Departures  from  cube-root  scaling  in  large  events  are  sometimes 
attributed  to  the  importance  of  gravity  forces.  For  surface 
bursts,  however,  departures  from  cube-root  scaling  have  also 
been  observed  with  small  charges  where  the  gravity  forces  are 
negligible  compared  to  the  material  strength.  The  apparent 
crater  data  seem  to  require  a value  for  n2  different  from  ni , 
which  would  violate  geometric  similarity. 


(Oh  recent  studies  in  the  Mine  Shaft  series 

t24,41)have  indicated  that  true  crater  dimensions  for  surface 
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explosions  in  rock  scale  fairly  consistently  for  charge  weights 
ranging  over  five  orders  of  magnitude,  but  not  by  cube-root^ 
scaling. ^ The  best  correlation  was  obtained  by  scaling  the  true 
crater  dimensions  by  height  of  burst  by  ^ ^ 

igure  15  shows  true  crater  dimensions  as  a function  of  charge 
weight  for  true  surface  shots  (dashed  curves)  and  HOB  = 0.9  r 
(s.olid  curves),  showing  remarkable  agreement  with  the 
W-*-'  ' scaling. 

n Experiment  Station  compendium  (Reference 

and  pp.  132  to  135)  presents  graphs  and  tables  in- 
uding  some  of  the  same  data  in  a different  fashion.  A best 
eight  rather  scattered  available  data  points  for  basalt 
and  granite  (based  on  two  charge  sizes)  for  scaled  HOB  between 
0.05  and  0.20  ft/lb^'^  was  represented  by  r-^-  = 0.195 
and  d^  = 0.023  wO-518  while  eight  points  (seven  charge  sizes) 
for  scaled  HOB  between  -0.05  and  +0.05  ft/lb^/^ 
presented  by  r^  = 1.211  wO-315  and  d^  = 0.588  wO-271^  Crater 

Donen+^?/9  7^-  for  W in  pounds  of  TNT.  Hence  the  ex- 

ponent 1/2.7  - 0.37  suggested  by  Figure  15  is  by  no  means  universal. 

^ For  the  purposes  of  the  present  study  the  two  calibration 

chaJgL""^^he  "he  la^geS 

Charges.  The  calibration  series  consisted  of  one-pound  spheri- 

charges  on  one-yard  cubical  blocks  of  cement  grout  and  1000- 

F?an^  spherical  charges  in  a natural  granitic  rock  medium? 

eiSnt  Reference  24  reporting  on  the  Mine  Ore 

two  100-ton  spherical  charges . The  figure  also  includes 
ata  from  two  comparable  series  of  previous  true-surface-shot 
cratering  tests.  The  Flat  Top  I event^^^^^  had  20-ton  TNT  spheres 

Multiple  Threat  Cratering  Experiment?*^^! ^ 
marked  MTCE  in  Figure  15,  had  4000-  and  16,000-pound  TNT  detona- 
tions in  basalt,  four  of  which  were  spherical  charges  in  geomet- 
ries comparable  to  those  of  Mine  Ore.  cnurges  in  geomet- 

19  one-pound  calibration  shots  were  especially  inter- 
esting, since  the  cement  grout  had  compressive  strength  and 
pulse  velocity  comparable  to  those  of  Lrmal  conSe?l  Tht 
average  static  compressive  strength  was  about  4500  psi.  Table 
3 summarizes  the  properties  of  the  grout.  This  serLs 

f?om^H0B^-^n  of  varying  the  height  of  burst 

u P • = 2rc  with  the  results  shown 

for  tills  h 
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TNT)  shots ‘-D  (,3,^  ^ 

to  -0.3^  ?or  ?ho=  ® f ”°®  '’®"S“S  ft.om  2rc 

was  beloS'the  rocrsur?I«  (HM  f 1 S 

hole  was  excavated  rrma>tit  -p  ^ charge  emplacement 

layer  ol  high!Slr?v-s?reittrj!fr”i"«  =»>^^8e  shape.  A 

the  hole.  shapjrwitfr?fSla?r?o  «^\Plastered  around 

to  harden.  Hydrostone  oemeS  wJ<-  ah's  allowed 

mate  coupling  from  tS^  Sge  ?^tte  nranitf  inti- 

sho?°"tte  "h"  ^ "-OPo®  “oJ'the  ^i??a:e-?f"  ent 

ze?o%GS)®^rin„=rd:?LlTi»ger"fh:n^?if 

the  joint.  Some  parts  of  thP  perpendicular  to 

planes  due  to  the  jointing  Desoite^th^’^  boundaries  were 
crater  dimensions  varied  fn  a coListSt^fa^S^  ’ average  true 
cept  for  one  shot.  The  combinenes^!?o  ^^7 

deptranfcJS?erSufon''S5?)^f^  th/dependencrorcJa??^'''' 

in  Figure  15  are  shown  in  a^^d^ 

The  upper  Figure  16  plots  true  crater  dan-v,  a • 
radii  versus  HOB  also  in  charee  r ^ charge 

equivalent  to  conventional  cube  root  scaling  ^'^li^th 

piiur:  i\®  sLieS  hy 

radius  rtw!/?Tf  ^rs^sloB/wf/t'^J^^'^^h?!  ' 

"Since  the  position  of  the  bottom  of'^a^qnh  "r^^^nted  as  follows: 
respect  to  the  medium  surface  i?  charge  with 

hob's,  the  HOB  must  stixl  be  scaled  h range  of 

the  relation  between  the  phvsical  qn^^-^  ^ ^rder  to  preserve 

the  bottom  of  the  chargS  tho  charge  center, 

addition,  the  curv^drSn  medium...  in 

acquires  an  'S'  shaped  bend\heS  the  HOR^  P^mts  (Figure  17) 
the  surface  tangent  geometrv  Th-;!^  2°^  ^ close  to 

the  cratering  potential  of  an  HE^nhir,°^  course,  indicates  that 
to  the  position  of  the  bottom  of  sensitive 

IZltty^  ^han  to  the  onlrittl  Tf  SI 

only  Ljrsur^lcHange^o?  hSe'^ 

also  suggest  the  need  for  studviL  1-L observations 
other  than  spherical  in  order  to  draw 


TABLE  3.  PROPERTIES  OF  CEMENT  GROUT 
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TABLE  4.  TRUE  CRATER  DIMENSIONS 


All  charges  were  1-pound  TNT  spheres  with  a radius  (r^,)  of  1.63  Inches, 


Charge 

Position^ 


True 

Crater 


True  Crater  Radius,  r^ 


.True 

Crater 


Shot 

No. 

HOB 

Tc 

dt 

fFeet) 

Maximum 

(Feet) 

Minimum 

(Feet) 

Average 
( Feet ) 

(In3) 

1 

0 

0.23 

0.74 

0.59 

0.68 

/'  240 

2 

0 

0.23 

0.74 

0.61 

0.6S 

' 200 

Avg 

0.23 

Avg 

0.74 

Avg 

0.60 

Avg 

0.67 

Avg  220 

18 

0.5 

0.27 

0.70 

0.55 

0.60 

_ _ 

Avg 

0.27 

Avg 

0.70 

Avg 

0.55 

Avg 

0.60 

3 

0.7 

0.17 

0.60 

0.45 

0.49 

no 

4 

0.7 

0.17 

■0.58 

0.37 

0.47^ 

70 

17 

0.7 

0.16 

0.52 

0.40 

0.45 

Avg 

0.17 

Avg 

0,57 

Avg 

0.41 

Avg 

0.47 

Avg  90 

5 

0.8 

0.14 

0.43 

0.30 

0.37' 

34 

6 

0.8 

0.17 

0.60 

0.39 

0.50 

16 

0.8 

0.11 

0.47 

0.33 

0.39 

Avg  0.14  Avg  0.50  Avg  0.34  Avg  0.42  Avg  34 


Avg 


0.09 

0.10 

0.10 

0.10 


Avg 


0.48 

0.43 

0.45 

0.45 


Avg 


0.26 

0.24 

0.21 

0.24 


Avg 


0.38 

0.31 

0.28 

0.32 


Avg 


22 

22 


0.07 

0.09 

0.08 


0.22 

0.22 

0.21 


0.19 

0.19 

0.13 


0.20 

0.21 

0.17 


7.8 


Avg 

0.08 

Avg 

0.22 

Avg 

0.17 

Avg 

0.19 

Avg  7.8 

Avg 

0.09 

0.09 

Avg 

0.68 

0.68 

Avg 

0.47 

0.47 

Avg 

0.57 

0.57 

Avg  — 

Avg 

0.08 

0.08 

Avg 

0.33 

0.33 

Avg 

0.13 

0.13 

Avg 

0.23 
0.23  ■ 

Avg  — 

Avg 

0.07 

0.07 

Avg 

0.32 

0.32 

Avg 

0.15 

0.15 

Avg 

0.24 

0.24 

Avg  — 

IN  CHARGE  RADII 
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Figure  17.  True  Crater  Radius,  Scaled  by  ^ Function 

of  Height  (or  Depth)  of  Burst,  Scaled  by  wl/3  (Geometric  Scaling), 
for  Spherical  TNT  Charges  Fired  in  Hard  Rock 


about  cratering  from  surface  bursts  of  HE.  Nevertheless  th^^ 
results  presented  establish  the  fact  that  cra?erlng 

sSrIacrcontac?  cSndi??onj!"  controlled 

for  the  1000-pound  charges  in  the  HOB  rLge  above  Ihe  surfacf 

shStrin'^this^Hnr''^"'*  craters  formed  from  the  1-pound 

shots  xn  this  HOB  range  appeared  to  be  spall  craters  resulting 

shock  stresses.  Two  possible  explanations 
eno^^of  were  suggested:  (1)  the  possible  exist- 

that  another  scaling  relation  for  these  HOB's,  or  (2) 

strength  (as  opposed  to  en  masse  strength)  of  the 
fhJ  so  much  higher  than  that  of  the  grout.  Although 

masse  strength  of  the  granite  is  usually  low  due  to^ 

^ ^ cracks , weathering,  etc.,  care  was  taken  in  the  selec- 
tion  of  the  calibration  series  test  areas  to  provide  competent 

fS^efbrairMLilir 

= 1 V ^^^“l^st-induced  stresses  in  the  1000-pound  shots  mav 

simply  be  a result  of  higher  medium  strength.  ^ 

ooncrete^ira''s!r?c,'^Jf''^®'*J‘*  ??  develop  empirical  formulas  for 
Tha  ^ ^ controlled  model  studies  all  in  concrete 

thL^tt^sri^thf  accurately 

tnan  those  in  the  field  tests  over  granite.  Small-scale  tests 

e.g.,  with  charge  masses  in  the  ratios  1:8:27,  could  be  used  ’ 

5"the1ca1e°ef?ec?f'“  and ’soSf  LsLsmSt 

on  extending  the  model  to  largL“chargfr?anrmore''Sost^y^^^^ 
ced^rttl;  in?erp??a?Ion?  estimation  pro- 

creta^Li^®  ‘T?  =“81  ested< ‘‘‘t  > that  the  complex  behavior  of  con- 
?eason  ?he  = "^self  to  scaling  procedures.  One  possiWe 
reason  for  the  suggestion  is  that  concrete  exhibits  rate  denenri 

elasticity  can  increase  brso  peJeen?  aS^ 

increaserftofS?iO-3lo"?oi^^c-J!"^N:ve“Lr""“-"^"? 

at  least^as^^^  scaling  in  soils  and  in  rock  media^that  are^ 

^ i complex  as  concrete  in  their  behavior,  there  is  a 

probability  of  reasonable  success  with  a moderate  effort. 

±n  following  section  will  survey  briefly  some  approaches 

the  prediction  of  dynamic  deformation,  penetration  and  ora 
concrete  that  are  of  a more  fuAdLenta?  cDac?er 
than  the  empirical  scaling  formulas. 
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3 . 3 More  Fundamental  Approaches 
3.3-1  Introduction 


blockf'^d^J^^hf  effects  on  finite  concrete 

2 ^ ^ (including  penetration  by  shaped  charges 

is^to^hP  explosive  cratering,  spall  and  corLr  fractures) 

i-t7r,ir.  "treated,  a complete  analysis  may  be  needed  of  some 
Th?^  order  to  establish  practical  guidelines 

This  section  sketches  the  state-of-tL-art  as  reJe^led  bO 

t leboratory  visits.  A great  deal 

has  already  been  accomplished  in  this  difficult  SL  but^ 

research  is  needed  before  it  will  be 
possible  to  prescribe  a specific  series  of  computL  simulJ 
tions  of  surface  and  near  surface  shots  (incJSdinrpenetral 
Which  would  be  reasonably  certain  to  furnish  a 
set  of  tables  or  curves  that  could  be  used  for  practinal 
damage  estimation  (something  comparable,  for  exLpl^to^he 
Lawrence  Livermore  Laboratorv's  At  He  fr.4.  ro4-* 
of  Underground  Nuclear  KvL  M°’^  Estimatinfi  Effects 


of  inreJtiS?JSnf interrelated  types 


1. 


Computational  (Code  development) 


2.  Theoretical  (Constitutive  equation  formulation) 


3.  Experimental  (Determination  of  properties  and 
verification  of  predictions). 


i::i?se\°sr?Se  x?i:rir:Sff  -s^r^seim 

?^r-d-  d™  e f 

since  they  are  quite  nerentTokL"!^  SciS®^”’ 


the  most  extensive  experimental  wav;  propagitjof datf ?hen 
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available  wasby  Gregson^Sl),  ^ determine 

uniaxial  strain  response  of  thin  plate  specimens  to  flyer 
plate  impacts  (stresses  from  40,000  psi  to  8 000  nnn  ncn  ■> 

a two-wave  structure  was  confirmed  and  a relatively  dispersive 
c,ne^t'^  Precursor  was  observed.  Read  and  MaidercLpa?ercaI! 
ulations  based  on  two  constitutive  models  (a  simple^elastic- 
plas tic  model  and  a porous  medium  model)  with  Gregson’s  ex- 

’ neither  gave  good  agreement.  They  concluded 
that  additional  experiments  are  needed  to  characterize  the 
Hugoniot  in  the  0 to  720,000  psi  range  which  it  o?thJ  Stat- 
t practical  importance.  Moreover,  since  there  was  a pomT,^a.^. 

attenuation  of  thin  sttett  ptlttt,  ?ht?  cot! 
attenuate Systematic  studies  of  release  wave  propagation  and 

There'^it  ttiu''i°?aS'’tf  1”1. 

u-_  . ® ^ of  the  kind  of  experimental  data  needed 

“"iaxial  strain  response  of  ootttet; 
to  Short  duration  pulses,  despite  the  great  importance  of  surh 
information  in  the  area  of  defense  analysis  and  deSgL 

uniaxial  strain  conditiciiS  prevail  ir.  V i3<a+ 
ioadtn®'  tnat  act  over  a large  surface  area  tnfaIJo  if  the 
loading  spherical  wave  front  at  large  distances  from  a noint 

expertise  for  the  uniexill 

fiwo-n  (eitnex  by  gas  gun  or  explosively  driven 

flyer  plates)  are  so  well  developed,  it  should  be  used  to  obtain 
he  needed  data  for  a number  of  well-defined  concrete  materials. 

For  application  to  near-surface  shot  explosive  loadincr 
penetration  problems  it  will  be  necessary  to  go  beyond  the 
untaxtal  strain  studies.  The  next  step  will  bf to  coSside? 

tSrf  if  if'"®?'  divergent  nature  of  thefaves 

pffp-u.  tendency  to  shock  up  and  consequently  greater 

Inf  f f Shear  strength  than  in  the  uniaLS  case 

tioffftfhP  f would  give  important  additional  informa- 
tf f penetration  behavior  and  alL  the  response” 

explosive  loading  over  a small  area.  These 

symmetry  but,  at  best,  axial 

searcrmJy  brSftlnguJsLd?"  analytical  re- 
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(1)  Uniaxial  strain  waves 

( 2 ) Spherical  waves 

(3)  Axially  symmetric  problems. 

A great  deal  of  work  has  been  done  in  the  first  configuration 
with  other  materials  but  very  little  with  concrete.  The 
spherically  symmetric  case  has  recently  been  the  subject  of 
some  studies  in  cement  grout,  primarily  for  testing  in  material 
gage  recording  techniques  to  be  used  later  in  rock  media.  The 
axially  symmetric  case  presents  additional  problems  both  experi- 
mentally and  computationally.  It  involves  two  space  dimensions, 
greatly  increasing  the  number  of  experimental  recording  positions 
needed  to  describe  the  wave,  and  it  admits  shear  cleavages  or 
sliding  surfaces  that  complicate  the  computation.  The  axially 
symmetric  geometry  will  eventually  have  to  be  treated  for  prac- 
tical applications , but  it  would  seem  more  profitable  at  this 

knowledge  to  pursue  the  spherical  case  first.  Section 
3.3-2  discusses  some  approaches  to  this  case  that  have  been  made 
or  proposed.  Then  Section  3.3-3  considers  constitutive  equations. 


3.3-2  Spherical  Wave  Experiments 

R.  P.  Swift ^52)  Qf  Physics  International  has  proposed  an 
experimental  procedure  for  spherical  wave  loading  of  large 
concrete  blocks  by  buried  charges,  using  in-material  piezo- 
resistive  stress  gages  and/or  electromagnetic  particle-velocitv 
gages  to  record  the  wave.  Figure  18  shows  one  of  the  several 
configurations  that  Swift  suggested.  By  placing  the  explosive 
sphere  off-center  in  the  block,  the  response  at  several  differ- 
j;adii_can  be  monitored.  Physics  International  has  applied 

Iif(53)  K ^ of  geologic  media,  including  gran- 

ite , tuff,  basalt,  clay,  and  sandstone^®^^ . Application  to 
concrete  would  be  straightforward  and  would  yield  important 
information  If  the  concrete  mix,  cure,  etc.,  are  well  defined 
and  controlled.  Figures  19  and  20  show  some  of  the  records  in 
granite  and  sandstone.  The  particle  velocity  gages  are  believed 
to  be  somewhat  better  than  the  stress  gages ,bSt  both  Ire  good 
enough  to  provide  useful  results.  ^ 

Stanford  Research  Institute  a preliminary  program  used 
concrete  blocks  for  testing  in-material  gage  measurement  tech- 
ques  for  later  use  in  soft  rock^^S).  Concrete  was  selected 
because  it  provided  an  inexpensive  and  convenient  experimental 
medium. ^ Figure  21  shows  the  test  geometry.  The  liquid  explo- 
sive  (nitromethane)  was  introduced  into  the  aluminum  sphere 
through  the  2-inch  access  tube  about  15  minutes  before  the  shot 
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11-A  of  Reference  55.) 
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The  particle  velocity  gages  were  of  a different  type  (mutual 
inductance  gages,  introduced  by  Engineering  Physics  Company ^5 6) 
and  adaptable  for  use  in  the  field).  Figure  21  shows  an  assembly 
with  three  stress  gages  and  three  velocity  gages.  Six  blocks 
were  tested  in  the  series  reported  in  Reference  56. 

Some  of  the  gages  were  actually  cast  into  the  material, 
while  others  were  inserted  into  holes  that  were  then  filled 
with  grout.  In  all  six  experiments,  concrete  cure  was  about 
two  weeks.  Drill  hole  gages  were  installed  seven  days  after 
the  pour,  and  the  grout  in  the  holes  cured  for  about  a week  be- 
fore the  shot.  Considerable  difference  was  noted  between  the 
in-material  gages  and  the  grouted  gages  in  the  measured  impulse, 
ihe  stress  from  the  in-material  gages  decayed  shortly  after  the 
peak  (e.g.,  20  to  30  microseconds)  but  remained  large  for  the 

200  microseconds.  There  was  some 
difference  in  the  grcut  mix  and  that  of  the  parent  concrete, 
and  possibly  the  surrounding  porous  concrete  drew  moisture  out 
of  the  grout  during  its  cure.  Possibly  a better  matching  grout 
can  be  found,  but  the  experiments  show  the  necessity  for  care- 
ful validation  of  any  proposed  in-material  gage  technique. 

With  metals  and  some  other  materials  in  hypervelocity  im- 
pacts a considerable  amount  of  information  is  obtainable  by 
measurements  of  back-face  velocities  and  stresses  by  various 
techniques  on  specimens  of  different  thicknesses.  Classic 
techniques  using  fly-off  pellets  of  the  same  material  as  the 
specimen  can  be  used,  and  should  in  particular  by  advantageous 
in  validating  the  results  from  in-material  gages. 

^ Prater^ has  used  a target  whose  basic  configuration 
IS  a half  cylinder  impacted  normally  at  the  center  of  the  flat 
face.  When  the  cylinder  radii  became  large  enough  so  that  peak 
stresses  fell  below  values  where  the  target  strength  could  be 
Ignored  (i.e.,  below  the  hydrodynamic  shock-wave  region),  the 
shape  of  the  target  was  changed  so  that  the  curved  surface  was 
replaced  by  several  flat  surfaces  each  perpendicular  to  the 
radius  drawn  from  the  impact  point  and  all  at  the  same  distance, 
tly-off  pellets  were  mounted  on  these  surfaces  as  well  as  piezo- 
electric  pressure  sensors.  A modified  larger  concrete  version 
of  such  a target,  loaded  by  a surface  explosive,  could  be  used 
to  measure  the  waves  emanating  from  the  crushed  and  cratered 
region  near  the  loaded  face.  This  gives  a situation  more  like 
the  surface  shots  considered  in  preceding  parts  of  this  section 
than  the  spherical  wave  geometries.  Apparently  no  such  experi- 
ments have  been  made  or  proposed  to  be  made  in  concrete,  but 
they  appear  to  be  feasible.  Preliminary  experiments  without  in- 
material  gages,  using  only  the  fly-off  pellets,  would  be  very 
easy  to  perform.  With  the  preliminary  results  available  to 
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guide  experimental  design  it  might  be  worthwhile  to  use  in- 
material  stress  and/or  velocity  gages  in  a subsequent  series 
of  tests. 

Section  3.3-3  considers  the  formulation  of  constitutive 
equations. 


3.3-3  Constitutive  Equations 

^ In  the  immediate  vicinity  of  the  explosively  loaded  area, 
loading  phase,  the  hydrostatic  pressure  is  so  large 
that  the  pressure-volume  equation  of  state  governs  and  material 
strength  par^eters  are  negligible.  The  study  of  many  other 
solids  in  this  hydrodynamic  regime  is  well  advanced.  For  sur- 
veys see  References  46,  58  and  59.  As  Read  and  Maiden(50) 

pointed  out,  however,  the  needed  data  on  concrete  is  still  not 
d.Vd.1  idble  • 

, unloading  phase,  even  for  the  regions  adequatelv 

described  by  the  hydrodynamic  theory  during  loading,  the  situa- 

Because  of  attenuation  by  overtaking 
spherical  divergence,  shock-wave  ampli- 
levels  where  material  strength  properties 
nn^+  dynamic  strength  properties  cannot  be  ade- 

quately measured  in  laboratory  experiments  with  homogeneous 
stress  states,  but  must  be  determined  in  situations  where  wave 
propaption  is  occurring.  As  a result,  these  properties  are 

largely  unknown  for  most  materials  and  certainly  unknown  for 
concrete « 

n + kinds  of  constitutive  equations  have  been  postu- 

lated, with  some  unknown  parameters  characterizing  the  material 
properties.  Predicted  wave  propagation  response  for  specific 
case.,  IS  then  compared  with  experimental  records.  The  computa- 
be  made  witn  assumed  values  of  the  parameters  and  then 
repeated  with  more-or-less  sophisticated  parameter  identification 

a best  fit  thS 

ShK  Iterative  procedure  can  become  very  expensive 

when  each  iteration  involves  a large  computer  solution  of  a wave- 

matprffli^°^  problem,  (tfg^fngian  analysis  based  on  multiple  in- 

veJsion^of  iterations.  One 

version  of  this  will  be  discussed  briefly  in  Section  3.3-4.  Some 

f the  constitutive  formulations  which  have  been  proposed  are 
discussed  in  the  following  paragraphs. 

A 1971  survey (64)  listed  four  inviscid  and  isothermal  mathe- 
matical models  for  geologic  materials , representing  the  then 
current  atate  of  the  art:  (1)  an  ideal  elastiSipJfsWo  S?Lial 
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model  with  different  constants  in  loading  and  unloading;  (2) 
a variable  moduli  material  with  different  material  constants 
in  loading  and  unloading;  (3)  a capped  elastic-plastic  material 
for  soils  in  which  the  cap  on  the  yield  condition  is  made  a 
function  of  the  plastic  volumetric  strain;  and  (4)  a capped 
elastic-plastic  model  for  rock  in  which  the  cap  is  a function 
of  a strain  parameter  (strain  hardening).  The  advanced  elastic- 
plastic  material  models  were  developed  from  the  Drucker- 
Prager  soil  mechanics  model,  with  deviatoric  yield  stress  a 
function  of  the  mean  stress.  This  led  to  a nonlinear  yield 
envelope  in  stress  space,  whose  normal  had  a component  such 
that  (by  the  associated  flow  rule)  a volume  increase  or 
dilatancy  always  accompanied  shear  deformation.  These  models 
can _ approximate  much  of  the  experimental  data  and  are  the 
easiest  to  fit,  but  they  cannot  match  standard  triaxial  tests, 
which  lead  to  compaction  instead  of  dilatation. 

The  variable  moduli  equations  resemble  differentiated 
isotropic^ elastic  equations;  S^j  = Gsij  for  the  deviators 

and  p = Ksj^j^  for  volume,  but  G and  K are  functions  of  the  dilata- 
tion (or  of  the  pressure)  and  G may  also  depend  on  the  deviatoric 
stress  invariants,  with  different  dependences  in  loading  and  un- 
loading. The  variable  moduli  models  are  relatively  easy  to  fit 
and  have  given  the  best  fit  to  most  of  the  data  available.  Be- 
cause they  do  not  incorporate  a yield  transition,  they  are  simple 
to  use  in  a computer  program.  As  in  deformation  theories  of 
plasticity  their  use  is  appropriate  only  for  proportional  or  near 
proportional  loading.  They  can  be  used  for  spherically  symmetric 
waves . - j j 

The  cap  models  represent  a further  development  of  the  ad- 
vanced elastic-plastic  models,  closing  the  yield  surface  on  the 
compression  side  with  a cap  whose  normal  direction  (over  most 
or  possibly  all  of  the  cap)  leads  to  compaction  instead  of  dila- 
tation. An  indirect  approach  is  required  to  fit  experimental 
data.  Most  applications  so  far  appear  to  have  assumed  an  ellip- 
tic shape  for  the  cap  in  the  ^ ~ plane.  is  the  first 

invariant  of  the  stress  and  is  the  second  invariant  of  the 
stress  deviator).  The  elliptic  cap  has  been  taken  with  center 
on  the  Jq  axis  and  with  horizontal  and  vertical  semiaxes  in  a 
constant  ratio.  Further  requirement  as  to  how  it  meets  the  yield 
envelope  (e.g.,  tangent)  then  leaves  only  one  parameter  to  fit 
to  hardening  data.  The  selection  of  the  elliptic  shape  for  the 
cap  IS  purely  arbitrary,  but  it  was  reasonably  simple  to  fit  and 
seemed  to  give  the  right  kind  of  shape.  A more  general  shape 
could  be  used,  depending  on  more  than  one  hardening  parameter. 

The  cap  models  appear  to  be  capable  of  matching  all  the  stand- 
ard test  data  and  much  of  the  wave  propagation  data. 
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The  models  discussed  above  are  all  rate  independent.  Ex- 
perimental evidence  suggests  that  material  response  is  rate 
dependent.  The  models  are  fitted  when  possible  from  experi- 
mental data  corresponding  to  loading  rates  comparable  to  those 
in  the  problem  of  interest.  Viscoplastic  models  can  also  be 
developed  when  enough  data  is  available.  The  most  primitive 
version  of  an  elastic-plastic  material  was  one  of  the  two  con- 
stitutive equations  used  by  Read  and  Maiden^^O)  to  try  to  model 
the  dynamic  uniaxial  strain  data  of  Gregson^^l)  for  concrete. 
This  model  assumed  no  hardening  and  no  plastic  volume  change 
and  gave  poor  agreement  with  the  experimental  results.  The 
other  model  used  was  a porous  medium  model,  which  also  gave 
poor  agreement.  S.  J.  Green,  et  al,  at  Terra  Tek,  Inc.  have 
recently  completed  a program'^S)  investigating  the  low  strain 
rcite  and  moderate  stress  level  of  behavior  of  concrete,  includ- 
ing a systematic  experimental  study  of  hydrostatic  and  multi- 
axial  stress  behavior.  More  advanced  elastic-plastic  models 
than  those  used  by  Read  and  Maiden,  including  the  cap  models, 
have  been  more  successful  in  rock  medial and  should  also  be 
suitable  for  concrete.  Although  it  appears  that  concrete  data 
have  not  yet  actually  been  fitted  to  any  of  these  advanced 
elastic-plastic  or  cap  models,  Swift(52")  in  1972  suggested  a 
cap  model  as  appropriate  for  modeling  the  data  he  proposed  to 
obtain  from  spherical  wave  experiments. 


It  is  known  that  rock  and  concrete  both  flow  when  subjected 
to  deviatoric  stresses  under  high  hydrostatic  pressures.  Thus, 
plastic  models  for  these  materials,  usually  considered  brittle, 
are  not  unreasonable  under  high  hydrostatic  pressures.  Near  the 
surface  of  the  target  and  in  the  later  stages  of  unloading,  where 
cracks  and  voids  begin  to  open  up,  these  models  might  not  seem 
so  appropriate.  Nevertheless  they  have  been  used  in  soils  and 
rock  media  with  some  success. 


A different  approach  to  modeling  both  ductile  and  brittle 


fracture  in  metals  under  hypervelocity  impact  has  been  used  by 
Stanford  Research  Institute^®6,67) , has  b< 


been  observed  that 
in  the  incipient  fracture  phase  of  the  deformation  (e.g.,  back- 
face  spall  of  a thin  plate)  where  voids  (in  ductile  materials) 
or  cracks  (in  brittle  material)  form  and  grow,  propagate  and 
coalesce,  the  usual  assumptions  about  metals  no  longer  hold. 

The  behavior  of  metals  in  this  regime  is  not  unlike  that  of 
granular  soils  or  fracturing  rock,  exhibiting  dilatancy,  for 
example,  and  other  types  of  behavior  not  usually  associated 
with  metals.  They  have  modeled  this  regime  in  metals,  both 
under  uniaxial  plate  slap  loading  leading  to  back-face  spall 
and  under  hypervelocity  impact  of  a steel  plate  by  a nylon 
sphere^®®^,  the  last  using  sophisticated  two-dimensional  computa- 
tions, giving  remarkable  agreement  with  experiments.  The  nylon 
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sphere  impact  led  to  cratering,  subcrater  fractures,  a phase- 
change  zone,  and  back  surface  fractures. 

The  modeling  of  the  nucleation  and  growth  of  the  voids  or 
cracks  as  rate  processes  is  described  in  References  66  and  67. 
Modifications  to  the  procedure  used  there  for  Armco  Iron  are 
described  in  Reference  69.  One  modification  was  to  assume  the 
nucleation  of  an  exponential  distribution  of  crack  sizes  at  each 
time  step.  The  form  assumed  for  the  distribution  was  considered 
reasonable  in  light  of  actual  measurements  in  rocks^' ^ . This 
modeling  technique  could  also  be  applied  to  penetration  and  ex- 
plosive loading  of  concrete.  It  appears  to  be  the  most  advanced 
modeling  technique  available.  A major  research  program  might 
lead  to  the  design  of  a set  of  computer  simulations,  with  ex- 
perimental input  of  data  and  some  verifications  of  predictions, 
that  would  permit  construction  of  a set  of  graphs  and  tables  for 
design  guidance. 


3.3-4  Lagrangian  Analysis 

In  the  introduction  three  interrelated  types  of  investiga- 
tion were  identified:  computational,  theoretical  and  experi- 

mental. The  interaction  of  the  three  is  most  evident  in  the 
relatively  new  approach  based  on  multiple  in-material  gage  re- 
cords'- ijDZ,b  ; which  is  known  as  Lagrangian  analysis.  This 
approach  was  originally  developed  for  plane  waves,  but  has  been 
extended  to  spherical  symmetry.  It  involves  multiple  gage  mea- 
surements of  both  stress  and  particle  velocity  at  stations  in- 
side the  medium,  and  subsequent  direct  deduction  of  the  operative 
constitutive  relation  without  iterative  computations.  The  version 
of  this  method  that  is  being  developed  at  Stanford  Research  In- 
stitute*- / ^ ; does  not  make  use  of  the  two  different  phase  veloc- 

ities (one  for  stress  and  one  tor  particle  velocity)  used  in 
previous  methods 

ili'^s'tt’ate  the  ideas,  consider  a spherically  symmetric 
wave  propagation.  The  Lagrangian  (initial)  and  Eulerian  radial 
distances  are  denoted  by  n and  r,  respectively.  Initial  and 
current  densities  are  pg  and  p.  Radial  particle  velocity  and 
stresses  are  and  ar»  while  (f,  = - Og  denotes  the  difference 

between  radial  and  tangential  stress.  (For  spherical  symmetry 
the  value  of  <j,  characterizes  the  dsviatoric  stress.)  The  con- 
servation laws  of  mass  and  momentum  and  the  definition  of  particle 
velocity  furnish  three  equations: 


6 6 


p ^ ^2  Uh 


(I)  • 


“r  = (n). 


(60) 


(61) 


(62) 


among  the  five  unknown  functions  (p,  u^,,  Uy,,  r and  <b)  of  fh^  -hwr, 

Scolplele^SlcaSL 

xt  lacKs  the  (unknown)  constitutive  equations. 

*5  multiple  in-material  gage  data,  portions  of  the 
sc^lutron  surfaces  or  = dr(ft.t)  and  u^  = UrO-lt)  are  oonsLJS?ed 
namerically  (Figure  22). 


The  known  equations  furnish 
r(/i,t)  - h + f^u^dt 


(63) 


e\/(hyt)  - £y(hftQ) 


- - 


(r>2u^) 


where  = i _ (p^/p)  is  volume  strain,  and 

•••5l!¥fl-Ssfe)  ] 


(64) 


(65) 


The  time  derivatives  needed  in  Equations  (64)  and  r>^-n  k 

shock-wave^frontp^  The^analysis  tt^n'give^valSef  of Je! 

strM 

constitutive  equations  are  then  fitted  to  this  data 
Note  that  iterated  computations  over  the  whole  field  arp  no+- 
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FLOW  SU ,^FACE 


^GAGE  PROFILES 


^Stress  or  Particle  Velocity  Histories  from  Multiple 
ages  are  Combined  to  Form  a Surface  in  the  Space-Time  Plane 


It  may  be  noted  from  Equation  (65)  that  in  experimental 
where  the  stress  and  particle  velocity  derivatives 
are  large  compared  to  their  difference  (e.g.,  in  the  rapid 

will  be  difficult  to  determine  (p  accuratelv 
The  flow  field  in  this  region  is  also  insensitive  to  ^h^Joic; 

equation  for  the  deviatoric  stress.  In  the 
nloading  regions,  where  the  derivatives  are  smaller,  the  flow 
IS  more  sensitive  to  the  deviatoric  stress,  and  thranalvsis 
provides  better  values  for  <p  where  <p  is  more  imporLn?!  ^ 


The  SRI  group  has  applied  this  procedure  to  laboratorv 
experimental  spherical  wave  data  obtained  by  Swift  o??Jysics 

Westerly  granltet53)  and  of  sLd- 
■;  with  some  success.  The  procedure  is  being  developed 
- rock  media  sponsored  by  the  Defense 

Nuclear  Agency ( 72 , 55 ) . For  such  in  situ  field  meLu?emen?s 
grouted  into  boreholes  in  large  expanses  of  rock* 
media,  the  procedure  will  be  subject  to  some  uncertainties. 


'A-.:  ' 


both  from  the  individual  measurement  difficulties  and  from 

;Sr?hr;i^T:  !!■=''  “"««ainty  „omd  be  ILTcLIIT 

with  the  application  of  the  method  to  studies  in  concrete 
under  controlled  conditions.  *=i.uaies  in  concrete 

impact  on  metal  plates(66,6f?  be  hypervelooity 

grangian  analysis  for  any  of  thSL  is  llmf^ia  • ‘‘®- 

diffioulty  of  obtaining  good  experiment^  in 

in  the  regions  of  intereit.  in-material  gage  data 


of‘‘conSI?e”BaJs^  Laboratory  Study  of  Dynamic  Tensile  Fracture 
3.4-1  Introduction 

fractu“y“Sc:rL"?hri::"^t“  “rner 

r^iL^e^rbi  ?L"d“:^rki 

design  dynamic  loadings,  the 

twice  the  static  uUiL?r?^S!?i  J selected  as 

tional  studies  related  to  strength ^ 3;,  Some  addi- 

concrete  havJ  beJn  ^poJ?ed  L properties  of 

those  studies  indica^e^hat  f^om 

Crete  increased  with  the  rate  of  ? ”*P^® s si ve  strength  of  con- 
the  static  value  at  an  inJeJlS-  about  twice 

Furthermore,  the  dynamic  elastic in/in/sec. 
imately  1-5  times  ttfs^a^lrvaLr  and 

of  the  material  was  found  to  inonM  impact  resistance 

Soma  additional  testr”  tl^e  s?a?lf  loading, 

crate  have  bean  reported  in  Raferenctr75"a“i6^"''?hese  tLllTs, 


conducted  on  impacted  laboratory-scaled  long  bar  specimens , 
of  varying  cure  time  and  properties,  have  shown  that  the  wave 
propagation  velocity  is  similar  to  that  of  igneous  rocks, 
that  the  dynamic  elastic  modulus  is  considerably  greater  than 
its  static  value,  and  that  pulses  of  short  duration  and  high 
rise  time  propagate  with  little  dispersion  and  some  attenuation 
depending  upon  the  indicated  strain  level.  Recent  studies^ 
indicate  that  a fracture  strain  energy  criterion  may  provide 
a meaningful  criterion  for  design  against  spall  and  corner 
fracture  phenomena.  In  addition,  some  recent  data  reported 
on  brittle  materials  other  than  concrete'"”^  indicate  that  a 
statistical  correlation  between  static  tensile  fracture  and 
compressive  pulse  fracture  in  long  bar  specimens  may  be  possible. 

At  the  back-face  of  a flat  slcib,  the  arriving  compressive 
wave  may  be  approximately  a plane  wave  of  uniaxial  strain, 
which  is  reflected  at  a free  back  surface  as  a tensile  wave, 
leading  to  tensile  fracture  under  conditions  of  approximately 
uniaxial  strain.  Corner  fractures  occur  under  more  compli- 
cated combined-stress  conditions,  but  in  the  low-stress  elas- 
tic regimes  brittle  fracture  is  often  assumed  to  be  governed 
by  the  maximum  tensile  stress  and  to  be  independent  of  the 
other  stresses.  The  range  of  stresses  and  rates  for  which 
this  maximxim-stress  failure  criterion  is  valid  has  not  been 
established.  In  the  range  for  which  it  is  valid,  strength 
determination  by  bar  tests , which  are  relatively  easy  to  per- 
form, should  be  useful. 

In  such  tests  an  impulsive  axial  load  is  applied  to  the 
end  of  a long  bar  by  a longitudinal  impact.  The  compressive 
wave  arriving  at  the  far  end  is  approximately  a plane  wave  of 
uniaxial  stress,  which  is  reflected  as  a tensile  wave.  This 
reflected  tensile  wave  is  superposed  on  the  oncoming  compres- 
sive pulse.  For  exactly  square  waves  and  perfect  reflection, 
the  superposed  waves  would  give  zero  stress  for  some  distance 
from  the  reflecting  end,  until  the  tail  of  the  oncoming  com- 
pressive wave  arrived,  when  a sharp  tensile  loading  would 
produce  fracture.  In  general  the  pulses  are  not  square,  and 
the  excess  of  the  tensile  wave  over  the  compressive  wave 
builds  up  gradually  until  it  exceeds  the  dynamic  tensile 
strength  at  some  section  of  the  bar. 

In  the  present  investigation,  information  on  the  thres- 
hold tensile  fracture  stress  of  concrete  bar  specimens  has 
been  sought.  Included  in  the  results  are  data  on  the  static 
tensile  and  compressive  strength,  and  the  influence  of  the 
geometrical  size  effects  on  impact  fracture  characterization. 


70 


b 

I 


3.4-2  Specimen  Preparation 


Long  bar  cylindrical  specimens  were  prepared  from  high- 
ear  ly-strength  cement  (Portland  III).  Specimen  composition 
and  static  properties  are  shown  in  TeO>les  5 and  6.  The  aggre- 
gates A-3  and  A-4,  supplied  by  NL  Industries  of  Edgar,  Florida, 
had  the  properties  shown  in  Table  6.  (^50  the  mean  diameter 

of  the  sieve  size  that  passes  50  percent  of  the  sample.  The 
uniformity  ratio  is  = Deo/Dio*) 

Concrete  was  mixed  in  quantities  sufficient  to  produce 
static  tensile  and  compressive  test  specimens  as  well  as  the 
long  dynamic  tensile  fracture  specimens.  Nominal  dimensions 
of  the  bar  specimens  used  in  the  test  program  were  0.75  and 
1.5  finches  in  diameter  by  18  inches  long.  Some  photographs  of 
typical  specimens  used  for  the  static  tension-compression  tests 
as  well  as  the  long-bar  specimens  used  in  the  dynamic  impact 
studies  are  shown  in  Figures  23  and  24. 

PVC  tubes  were  used  as  molds  for  fabricating  the  cylin- 
drical bar  specimens.  Spraying  the  interior  surface  with  a 
mold  release  agent  insured  a smooth  surface.  For  casting  the 
specimens,  the  PVC  tubes  were  held  in  a vertical  position  in 
a special  jig  fixture  with  the  lower  end  sealed  with  a rubber 
stopper.  The  concrete  mixture  was  added  in  2-inch  layers  and 
consolidated  by  tamping.  A second  stopper  was  then  used  to 
seal  the  top  end  of  the  tube  and  the  sides  of  the  split  PVC 
tubing  sealed  with  waterproof  tape.  During  the  initial  set- 
ting  period,  the  tubes  were  placed  in  a horizontal  position  emd 
allowed  to  set  for  three  days  in  a moist  cure  environment.  The 
samples  were  then  moist-cured  for  an  additional  period  of  four 
days  at  which  time  the  specimens  were  removed  from  the  holds 
and  axr  dried  in  a low  humidity  room* 


3.4-3.  Test  Procedure  and  Results 

Static  tensile  and  compressive  tests  were  run  on  the  prin- 
cipal  material  types  with  test  data  reported  in  Table  6.  The 
tensile  test  specimens  used  had  a nominal  diameter  of  0.75  in. 
by  3 in.  long.  Steel  disks  0.25  in.  thick  were  epoxied  to  the 
top  and  bottom  surfaces  of  the  specimen  and  eye  bolts  inserted 
ii;  r attachment  to  the  loading  machine  as  shown 

in  figure  23a.  For  the  compression  snecimens,  nominally  0.75 
in.  diameter  were  0.5  in.  long  Teflon® sheets  were  used  at  the 
loading  surfaces  to  minimize  friction  effects. 


71 


SPECIMEN  COMPOSITION 
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Figure  24.  Dynamic  Impact  Bar  Specimens 
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The  tensile  fracture  studies  were  performed  using  an  air 
gun  assembly  described  in  Reference  78  and  shown  in  Figure  25. 
Some  typical  dynamically  fractured  tensile  specimens  are  shown 
in  Figure  26,  and  typical  strain  profiles  for  impacted  concrete 
bars  are  shown  in  Figure  27.  For  the  present  tests,  1.0  in. 
and  2.0  in.  long  cylindrical  steel  impactors  0.485  in.  in 
diameter  with  hemispherical  noses  were  fired  at  capped  and  un- 
capped concrete  bars  with  varying  velocities  so  as  to  bracket 
the  threshold  tensile  fracture  stress. 

The  approximate  values  of  the  threshold  fracture  stresses 
were  estimated  by  the  following  procedure.  After  the  tests  had 
been  made  on  uninstrumented  bars  of  various  sizes  with  differ- 
ent aggregate  sizes,  several  strain-gage  instrumented  bars  were 
tested.  The  dynamic  elastic  modulus  E = pC2  was  determined 
from  measurements  of  wave  speed  C and  the  known  density  p.  The 
instrumented  bars  were  impacted  first  at  very  low  velocities  to 
measure  wave  speeds  and  then  to  higher  velocities  approaching 
the  previously  determined  threshold  impact  velocity  for  tensile 
fracture.  From  the  maxxmum  straxn  recorded  at  the  gage  nearest 
the  first  observed  fracture  location  the  incipient  fracture 
strain  was  estimated.  The  stress  was  then  calculated  by  Hooke's 
Law  with  this  strain  and  the  previously  determined  dynamic 
modulus.  The  incipient  tensile  fracture  stresses  estimated  by 
this  procedure  were  on  the  order  of  1.3  to  1.5  times  the  static 
tensile  strengths  given  in  Table  6.  This  is  smaller  than  the 
factor  of  two  customarily  assumed (7 3)  -^o  relate  the  dynamic  ten- 
sile strength  of  reinforced  concrete  to  the  static  ultimate  ten- 
sile strength. 
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a)  Strain  pulse  propagating  in  a 1.50  inch 
diameter  concrete  bar,  Bar  9,  1",  195  in/sec 
C,  A-3,  7,  20jj  sec/cm,  0.5  m volt/cm 


b)  Strain  pulse  propagating  in  a 1.50  inch 
diameter  concrete  bar.  Bar  9,  2",  167  iii/sec 

C,  A-3,  7,  20y  sec/cm,  0.5m  volt/cm 


Figure  27.  Typical  Strain  Profiles  for 
Impacted  Concrete  Bars 
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Conclusion 


Net  enough  actual  data  on  cratering  in  concrete  was  found 
to  permit  the  deduction  of  any  specific  formulas  or  graphs  for 
concrete  similar  to  those  presented  in  Section  3.2  for  certain 
rock  media.  Because  of  the  similarities  between  the  material 
properties  of  concrete  and  of  some  of  the  rock  media,  it  is  be- 
lieved that  the  same  methods  that  have  been  successful  in  rock 
media  could  also  be  used  for  concrete  if  enough  experimental 
data  can  be  obtained  under  controlled  conditions. 


Empirical  scaling  methods  similar  to  those  represented  in 
Figures  15  to  17  for  spherical  charge  true  surface  shots  (i.e. , 
half  buried)  or  near  surface  shots,  with  charges  in  intimate 
contact  with  the  concrete,  could  certainly  be  applied.  Because 
controlled  tests  on  concrete  would  have  fewer  randc...  variations 
in  the  medium  than  those  in  geologic  materials,  the  irode.Ling 
should  be  even  better  in  the  lower  part  of  the  charge  weight 
range  than  the  very  good  agreement  shown  by  the  Mine  Shaft  cal- 
ibration series  of  Figures  15  to  17. 


Small  scale  tests  with  charge  weights  varying  over  a range 
or  the  order  of  one  pound  to  27  pounds  could  be  used  to  deter- 
mine parameters  in  such  empirical  models.  The  one-pound  charges 
in  cement  grout,  -eported  in  Table  4 of  Section  3.2,  were  placed 
on  cubical  target  blocks  one  yard  on  a side. 


4.,  crater  formed  had  a radius  about  one-fourth 

of  the  block  width  and  a depth  about  a third  of  the  radius. 

This  suggests  that  blocks  of  smaller  depth  than  width  could 
be  used  in  the  program  if  they  were  set  in  a dense  soil  or 

minimize  backface  reflection  and  scabbing.  This 
should  be  verified  by  comparison  shots  with  the  smaller  charges 
Geometric  scaling  would  require  a block  three  yards  on  a side 
for  the  27-pound  charges.^  A follow-up  program  should  test  the 
empirical  formulas  determined  by  the  small-scale  tests  with  a 
limired  number  of  larger  charges. 


Effects  of  charge  shape  and  contact  conditions  woulc.  also 

tion  recognized  at  the  beginning  of  Sec- 

tion 3.2  that  successful  empirical  representation  of  controlled- 
condition  surface-burst  loadings,  will,  however,  be  of  limited 
value  for  practical  prediction  of  bomb  damage  after  partial 
penetration.  The  importance  of  the  contact  conditions  and  the 
-2urst  cannot  be  overemphasized.  Reinforcement  and 
ffuite  taiget  thickness  will  also  be  important. 
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methods  using  numerical  computation  of  a 
^Ictc  fcve.^t  (inducing  penetration  by  a shaped  charge  follow- 
y explosive  cratering,  spall  and  corner  fractures)  are 
certainly  possible.  But  the  development  of  constitutive  models 
for  use  in  such  computatinnp  is  ^ fu,  ti  t ? 

search  phase  This  fundamental  resejr^^  has  madi^romrsinrad- 
vances  in  other  materials,  including  geologic  materials  with 

fSnLr“r"coJcrete?  =«ort  has  been 


Much  more  funnoinental  research,  especially  experimental  and 
heoretical  work  on  material  characterization  (Section  33)  is 
needed  before  a specific  series  of  computer  siLSions  d typi- 
cal cases  to  obtain  sets  of  tables  and  graphs  for  practical  ^ 
damage  estimation  can  be  prescribed.  Such  useful  design  aids 
should  be  developed  when  material  characterizations 
ufficiently  good  to  use  in  the  computer  simulations  are  avail- 


experimental  lines  of  investigation  are  recommended  to 

an^roLher'^^Thr^''''^^^  properties  for  these  more  fundamental 
pproaches.  These  are  uniaxial  plane  strain  studies  bv  date- 

describerL  SectiM  Ta  studies  similar  to  those 

co?d-  Wiih  ehff  “=^"8  multiple  in-material  gage  re- 

core...  With  these  experimental  records  available,  one  of  the 
i constitutive  mnadcc  3 -•  • 


bf  fitted'^t^th  r discussed  in  Section  3.3-3  codd 

the  spherical  wave  data  (for  example,  a capped  elas- 

LctiorrS-r'^^After  the  Lagrangian  analysis  method  of 

itCTion  3.3  4.  After  successful  modeling  of  the  spherical 

problems,  the  more  difficult  axisymmetric 
penetration  an\  cratering  problems  could  be  attempted 
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APPUWDiX 

LITERATURE  SEARCH,  LABORATORY  VISITS,  TECHNICAL  MEETINGS 

In  conjunction  with  the  specific  areas  of  investigation 
noted  in  the  introductory  remarks,  several  literature  search 
strategies  and  research  information  sources  were  used  for  per- 
tinent unclassified  subject  references.  The  table  below  identi- 
ties the  source  agency  and  search  title  used  for  this  purpose. 

TABLE  A-1.  LIST  OF  LITERATURE  SEARCHES 


Agency 

North  Carolina  Science 
and  Technology  Research 
Center 

North  Carolina  Science 
and  Technology  Research 
Center  

Defense  Documentation 
Center  - Literature  Surve'i 

Defense  Documentation 
Center  - Work  Unit 
Summary  Current  Programs 

Smithsonian  Science 
Information  Exchange ^ Inc. 


Search  Title 

Effects  of  Pressure  Blast  Loading 
in  Brittle  Structural  Materials 


Dynamic  Behavior  of  Concrete 


Blast  Response  of  Concrete  and 
Brittle  Materials 

Dynamic  Behavior  of  Concrete 


Dynamic  Properties  of  Concrete 


In  addition  to  the  above  literature  searches,  trips  to 

te^sile^fil^^?^  research  dealing  with  material  removal  and 
fracture  of  concrete  have  been  made.  Also,  several 

?Sck  Snd  soil^meS?rh"®^^i''®  dynamic  effects  in  concrete, 

attended  and  contacts  with 
active  investigators  in  these  areas  pursued. 

ed  is^nclSS  below?''®^  visited  and  pertinent  meetings  attend- 


TABLE  A-2.  PLACES  AND  LABORATORIES  VISITED 
AND  CONFERENCES  ATTENDED 


Research  Laboratories 


Conferences 


Waterways  Experiment  Station 


Accoustical  Society  of  America 
Meeting,  29  Oct-2  Nov,  1973 
(Wave  Propagation  Section) 


Lawrence  Livermore  Lab 


Mechanisms  of  Explosion  and 
Blast  Waves,  12-16  Nov.  1973 


U.  S.  Army  Corps  of  Engineers 
(Livermore ) 


American  Society  of  Mechanical 
Engineers,  11-15  Nov,  1973 
(Rock  mechanics) 


Stanford  Research  Institute 


Shock  and  Vibration  Symposium 
3-7  Dec,  1973  (Impact  and  Shock 
loading  section) 


Physics  International 


Seventh  U.S.  National  Congress 
of  Applied  Mechanics,  3-9  June 


The  computer  searches  produced  little  information  directly 
applicable  to  the  cratering/material  removal  problem.  They  did 
produce  abstracts  of  papers  in  related  areas  as  listed  below  in 
three  categories. 


TABLE  A- 3.  CATEGORIES  OF  ABSTRACTS  FROM  LITERATURE  SEARCH 


Number  of  Abstracts 


Blast  Loading 


Fracture 


Material  Pronerties 


i 


r 

I 


Many  of  the  abstracts  in  the  blast  loading  category  pertained  1 

to  response  characterization  of  concrete  structural  configura-  I 

tions  such  as  arches,  frames,  slabs,  and  beam  elements,  while  | 

those  categorized  under  fracture  were  associated  with  the  topic  t 

° fracture  of  materials,  including  concrete.  The  ab-  J 

stracts  listed  under  material  properties  were  concerned  mainly 
with  uniaxial  and  multiaxial  testing,  including  the  effects  on  t 

material  strength  of  such  variables  as  aging,  aggregate  particle  i 

size,  and  water-cement  ratio. 

References  cited  in  the  text  are  listed  in  the  following 
subsection.  Most  of  these  were  not  located  by  the  computer 
searches,  but  were  suggested  by  personal  contacts  with  active 
investigators  or  discovered  by  scanning  library  copies  of  cer-  * 

tain  journals.  • 
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